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I INTRODUCTION 

A.  The  Low-Angle  Tracking  Problem 

Radar  targets  at  long  range  and  low  altitude  are  very  difficult  to 
detect  and  are  especially  difficult  to  track.  This  is  because  the  radia- 
tion scattered  by  the  target  is  reflected  by  the  ground  to  the  radar 
antenna  and,  as  a result,  in  the  simplest  case  the  direct  return  is  mixed 
with  a coherent  scattered  return  (see  Figure  1).  This  scattered  signal 
adds  constructively  or  destructively,  depending  on  the  exact  geometrical 
relationships  between  the  target,  ground,  and  radar  antenna.  When  this 
process  occurs  at  very  low  radar  beam  elevation  angles,  the  scattered 
signal  falls  within  the  main  antenna  beam.  In  effect,  the  scattered  sig- 
nal is  a multipath  return  from  the  target  that  causes  an  "image"  of  the 
target  to  appear  below  the  horizon.  Thus  many  monopulse  radar  systems 
will,  when  tracking  a low-angle  target,  center  on  the  horizon  or  become 
somewhat  unstable,  depending  on  the  exact  phase  relationship  of  the  multi- 
path  and  direct  radar  return  signals.  At  times,  monopulse  radars  will 
even  track  the  subterranean  image  of  the  target. 

Techniques  to  improve  low-angle  tracking  capability  have  been  the 
subject  of  increasing  interest  as  general  radar  performance  has  improved. 

It  is  increasingly  important  to  know  the  limits  of  performance  possible 
in  ground-control ied-approach  and  enroute  air  surveillance  radar  systems 
as  the  air  traffic  around  airports  increases  and  the  danger  of  air  colli- 
sions increases. 

borne  very  expensive  weapon  systems  are  built  around  certain  assump- 
tions concerning  the  low-angle  tracking  ability  of  installed  defensive 
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rad.,,-  systems,  The  Payoffs  tor  gaining  . bettor  undorst„ndi„g  of 
angle  tracking  arc  further  examined  by  Smith  and  Moiling  (1974).  Barlon 
(1974)  has  recently  reviewed  the  lo-angle  tracking  problem. 


TARGET 


figure  1 


ILLUSTRATION  OF  THE  MULTIPATH 
FROM  A FLAT  EARTH 


PROBLEM  FOR  SPECULAR  REFLECTION 


u. 


electromagnetic  Wave  Scattering 

A key  issue  in  lew-angle  multipath  propagation  is  the  scattering 
process.  A,  very  low  grazing  angles,  the  scattering  process  is  much  more 
mplex  than  at  near  normal  incidence.  Shadowing,  multiple  reflections, 
suriace  roughness,  and  earth  curvature  effects  all  become  important  con- 
siderations  at  grazing  incidence. 


Because  of  these 
applicable  to  these  lc 


difficulties,  very  little  scattering  theory  is 
W angles'  Meaningful  experiments  in  these  low- 
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angle  regimes  are  also  very  difficult  to  perform.  The  result  is  that 
there  is  only  a little  information  available  that  concerns  the  single 
most  important  feature  of  the  low-angle  tracking  problem.  Section  1 1 of 
this  report  briefly  reviews  the  present  state  of  knowledge  of  low-angle 
scattering  processes  and  compares  some  relevant  lunar  bistatic  radar  data 
with  theory. 

c>  Radar  Systems  for  Low-Angle  Tracking 

Recently  several  radar  techniques  have  been  proposed  to  improve  low- 
angle  tracking  capability.  The  most  impressive  work  demonstrated  to  date 
is  that  of  White  (1974).  His  method  involves  complex  signal  processing 
of  the  radar  return  as  received  from  several  horns  of  a dish  antenna  or 
in  conjunction  with  an  array  antenna.  Dax  (1973)  proposed  a similar  ap- 
proach. Sherman  (1966,  1971),  Peebles  and  Berkowitz  (1968),  Peebles  (1971), 
Peebles  and  Goldman  (1971),  Howard  et  al.  (1971,  1973),. and  Symonds  and 
Smith  (1973)  all  proposed  and  reported  on  the  use  of  complex-angle  signal 
processing  of  the  radar  returns  for  low-angle  and  multiple-target  radar 
systems.  Other  signal-processing  approaches  have  been  suggested  by 
Von  Schlachta  (1973),  Sklar  and  Schweppe  (1964),  Pollon  (1967),  Pollon 
and  Lank  (1968),  and  Ksienski  and  McGhee  (1968),  These  are  more  general 
clutter  and  mul tiple- target  studies  but  are  applicable  to  the  low-angle 
trackit. g problems.  These  and  some  additional  possibilities  for  improving 
low-angle  tracking  capability  are  considered  in  Section  III, 

Section  III  also  considers  several  possibilities  for  rejecting  the 
reflected  radar  return  by  the  use  of  advanced  antenna  techniques,  terrain 
modification,  selective  siting,  and  radar  fences. 

Conclusions  and  recommendations  are  presented  in  Section  IV. 
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II  ELECTROMAGNETIC  WAVE  SCATTERING— APPLICATIONS 
TO  LOW-ANGLE  TRACKING 

A . Genera  1 

The  analysis  of  electromagnetic  wave  scattering  from  rough  surfaces 
has  become  increasingly  important  with  the  advance  of  radar  technology 
since  W'orld  War  II.  Theoretical  work  has  tended  to  bo  separated  into 
two  groups:  (1)  the  traditional  radar  community  with  interests  in  pre- 
dicting the  statistical  properties  of  the  scattered  signal  when  the  sur- 
face properties  are  fairly  well  known  (e.g.,  the  radar  clutter  problem), 
and  (2)  the  astronomical  and  geophysical  community  where  one  wants  to 
obtain  statistical  information  on  an  unknown  surface  from  a knowledge  of 
the  scattered  field  (e.g.,  radar  astronomy  of  planetary  surfaces  and  re- 
mote sensing  of  the  terrestrial  environment).  Information  from  both  these 
communities  should  be  sought  in  attacking  a particular  problem.  Some  use- 
ful general  and  review  sources  are  Barrick  and  Peake  (1967),  Barrick  (1970), 
Bass  and  Fuchs  (1972),  Beckmann  and  Spizzichino  (1963),  Evans  (1970),  Evans 
and  Hagfors  (1968),  Simpson  (1973),  Skolnik  (1970,  Chapters  22,  25,  and  26), 
and  Tyler  and  Ingalls  (1971). 

This  section  will  review  some  necessary  general  scattering  theory, 
but  will  concentrate  on  those  features  relevant  to  the  problem  of  low- 
angle  radar  tracking.  The  magnitude  and  statistical  nature  of  the  scattered 
component  at  the  radar  will  be  the  main  item  of  interest.  Several  questions 
arise.  For  example,  will  the  scattered  component  phase  be  coherent  or  in- 
coherent with  respect  to  the  direct  return,  and  on  what  time  scale  will 
any  phase  coherence  persist?  How  much  will  the  scattered  return  be  Doppler- 
shifted  with  respect  to  the  direct  return?  The  treatment  below  will  certainly 
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be  far  from  complete.  Our  aim  is  simply  to  report  the  features  of  scatter- 
ing theory  that  seemed  most  relevant  during  the  1974  JASON  summer  study. 

9.  Rough-Surface  Scattering  Theory 

1 . Brief  Overview  of  Scattering  from  Rough  Surfaces 

The  basic  scattering  geometry  is  illustrated  in  Figure  2.  A 
signal  is  transmitted  along  r toward  a target,  T.  It  is  scattered  by  the 
target  and  returns  to  the  radar  receiver  via  the  direct  path,  TR,  and  via 
a number  of  scatter  paths  such  as  TOR.  The  resultant  E field  at  R is  then 
the  vector  sum  of  the  wave  E fields  arriving  via  the  various  return  paths. 


FIGURE  2 SCATTERING  GEOMETRY.  A radar  pulse  transmitted  from  the  radar  site,  R,  is 
scattered  by  the  target,  T,  and  returns  to  the  radar  via  the  direct  path,  TR,  and 
a collection  of  paths  involving  a scattering  from  the  Earth's  surface  such  as  the 
path,  TOR.  The  heights  off  the  mean  surface  of  the  radar  and  target  are  hr  and 
ht,  ds  is  an  elemental  scattering  area,  and  v is  the  target  velocity  vector. 


Consider  a radar  receiver  that  transmits  and  receives  a horizon- 
tally polarized  wave.  As  a simplification  we  will  consider  only  the  single 

component  of  the  electric  field  parallel  to  the  surface  S and  perpendicular 
— » 

to  r.  This  is  not  unreasonable,  since  waves  arriving  at  R will  be  confined 
by  the  antenna  to  directions  very  nearly  along  r.  Other  electric-field 
components  can  of  course  be  considered  separately.  Let  the  E field  arriv- 
ing at  the  target  be 


e 


E o 
TR 


ikr 


The  diroct-ochc  E field  then  becomes 


ik2r  ifT 
E„e  p e 
TR  T 


RD 


(1) 


if  t 

where  P^e  is  the  reflection  coefficient  of  the  target  and  the  factor 

iwt 

e has  been  supressed  to  simplify  the  mathematical  expressions. 

To  calculate  the  electric  field  E arriving  at  the  receiver  via 

RS 

surface  scatter  paths,  we  will  use  the  Helmholtz  integral  of  physical 
optics.  At  each  pointj  0,  of  the  scattering  surface,  S.  we  will  calculate 
the  scattered  electric  field  by  assuming  the  existence  of  a reflecting 
plane  tangent  to  the  surface  at  the  point  0.  Following  an  approach  similar 
to  the  work  of  Beckmann  and  Spizzichino  (1963,  pp.  17-22  and  178-181)  we 
have 

E -1Uv1*T)eltr  ft  ylk(ri~2> 

ERS  2Ar  J J rir2 

S 


• (i+R) 


— » 
r 

1 


r 

1 


+ 


(1-R) 


ds 


(2) 


if 

where  R = p^e  s is  the  Fresnel  reflection  coefficient  for  the  reflecting 
plane  at  0,  n is  a unit  vector  normal  to  the  local  surface  at  point  0,  and 
gR  is  the  voltage  gain  of  the  receiving  antenna  along  r , This  particular 
formulation  of  electromagnetic  wave  scattering  from  a rough  surface  is 
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Since  the  surface  ma;  be  rough 


usually  called  quasispecular  scattering. 

A 

n is  not  necessarily  perpendicular  to  the  mean  surface.  In  fact,  the 
roughness  of  the  surface  is  modeled  by  letting  the  direction  of  the  local 
normal,  n,  vary  in  a random  manner  about  a mean  value. 


Though  there  are  alternative  formulations,  Eq.  (2)  is  especially 
useful  since  it  clearly  shows  the  actual  physical  summation  process  as 
well  as  other  important  general  features.  If  the  surface  S is  relatively 

homogeneous  and  smooth,  n will  be  very  nearly  perpendicular  to  the  mean 

„ ik(r1+r9) 

surface  and  the  exponential  term  e 1 * will  tend  to  dominate  the 

behavior  of  the  integral.  The  received  wave  will  thus  be  mainly  coherent 
and  similar  to  that  expected  from  the  specular-reflection  case  described 
below.  If  the  surface  S is  comparatively  rough  and/or  highly  nonhomo- 
geneous,  varying  randomly  from  one  scattering  point  to  another,  the  ex- 
ponential term  will  not  be  important  and  the  received  wave  will  tend  to 
be  incoherent.  Waves  scattered  from  rough  surfaces  in  nature  will  in 
general  have  both  coherent  and  incoherent  components.  It  is  the  coherent 
component  that  is  the  most  devastating  to  radars  tracking  low-elevation 
targets,  because  it  can  add  in  antiphase  to  the  direct  radar  echo  and 
cancel  the  desired  signal  completely. 


When  the  surface  is  homogeneous,  but  still  comparatively  rough, 
the  geometric  factor  (l/r^rg)  comes  into  play.  This  factor  tends  to  divide 
any  incoherently  scattered  return  into  foreground  (r^,  small)  and  horizon 
(r2>  small)  components  over  regions  where  g^  is  not  a strong  function  of 
the  location  of  the  scattering  point.  Barton  (1974,  p.  691)  and  Beckmann 
and  Spizzichino  (1963,  Chapter  2)  give  examples  of  specific  cases. 

Scattering  from  a rough  surface  is  best  understood  by  first  con- 
sidering a perfectly  smooth,  homogeneous,  infinite  flat  surface.  Reflection 

from  such  a surface  is  governed  by  Snell's  law  (0  = 6 , cp  = 0 in  Figure  3) 

i s s 
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and  the  classical  Fresnel  reflection  coeff ieients  given  in  Eqs.  (22)  and 
(23)  below.  Such  a reflection  (6.  = 0 , cp  = 0)  is  called  specular  and 
the  term  reflection  is  used  since  no  power  is  scattered  in  other  directions. 
Specular  reflection  is  cohe rent  in  that  the  reflected  field  parameters 
(amplitude  and  phase)  are  uniquely  predictable  — i.e.  , deterministic,  as 
shown  below.  Any  real  surface  is,  of  course,  finite,  so  some  energy  is 
scattered  into  other  directions  close  to  the  specular  direction 

(9s  = ®i’  cps  = but  for  thc  large  reflecting  areas  considered  here  the 
specular  component  will  dominate  strongly. 


FIGURE  3 LOCAL  FRAME  OF  REFERENCE  FOR  AN  APPROXIMATELY 
FLAT  SCATTERING  SURFACE.  The  angle  of  incidence,  0, 
grazing  angle,  y,  and  wave  vector  (k,  |kj  = 2nf\)  for  the  incident 
wave  are  denoted  by  the  subscript  i,  while  the  subscript  s identifies 
the  scattered-wave  parameters. 

In  considering  a rough  surface,  we  have  several  types  of  surface 
models  from  which  to  choose:  semi-empirical  models,  which  are  based  on 

simple  physical  ideas  and  matched  to  given  sets  of  data  by  parameter  ad- 
justments; geometric  models  in  which  simple  geometric  shapes  are  arranged 
randomly  on  a plane;  and  statistical  models  in  which  the  surface  height 
above  the  reference  (x,y)  plane  is  a random  variable.  The  statistical 
approach  is  the  most  general  and  we  shall  consider  such  models  for  the 
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most  part.  However,  a number  of  problcms--e.g. , diffuse  scattering-- 
have  not  been  satisfactorily  solved  via  the  statistical  approach  and  one 
must  use  an  alternative  model,  Barrick  (1970)  reviews  results  from  all 
three  approaches  with  an  emphasis  on  backscatter. 

In  the  statistical  approach  the  surface  is  characterized  by  an 
average  (rrns)  height  deviation,  h,  from  the  x,y  plane  and  a correlation 
length  £.  As  h increases  from  zero,  less  power  is  reflected  specularly 
and  more  power  is  scattered  into  nonspecular  directions.  A useful  param- 
eter is  defined  by  g = kh  (sin  y + sin  y ),  where  k = 2nA . For  a 

1 s 

slightly  rough  surface  (g  « 1),  the  coherent,  specular  component  will 
still  be  present;  but  an  incoherent  component  will  become  increasingly 
important  as  g increases, especially  in  directions  away  from  the  specular 
direction.  For  g ^ 1 the  surface  is  rough  enough  that  the  incoherent 
component  dominates.  Figure  4 illustrates  the  transition  from  smooth  to 
rough.  See  Beckmann  and  Spizzichino  (1963,  Chapter  5)  for  more  details. 


(a)  (b)  Ic)  (dl 

FIGURE  4 TRANSITION  FROM  SPECULAR  REFLECTION  TO  DIFFUSE 


SCATTERING.  The  surfaces  vary  from  a smooth  surface  with 
g = 0 in  (a)  to  a very  rough  surface  with  g » 1 in  (d) 

The  parameter  g is  defined  as  (27th/X)  (sin  y;  + sin  ys),  where 
h is  the  rms  surface  height  deviation  and  X is  the  wavelength 
(after  Beckmann  and  Spizzichino,  1963,  p.  90). 
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2. 


Specular  Reflection 


We  define  the  field  (as  opposed  to  power)  reflection  coefficient 

as  Beckmann  and  Spizzichino  (1963,  p.  22)  define  it:  p = E /E  where 

R/  RQ,  wneie 

ER  1S  thG  reCGlVed  flGld  and  Ero  the  field  that  would  be  reflected 
specularly  (Y  = Y^  = Yg)  by  a smooth,  perfectly  conducting  plane  under  the 
same  circumstances.  For  a perfectly  smooth  dielectric  surface,  p = p eJOr° 


where  p e 


3ao  t 


s given  by  Ramo  et  al.  (1966,  p.  358).  For  a slightly  rough 


surface  we  let  Pq  - P spec  [cf.  Beckmann  and  Spizzichino  (1963,  p.  80- 

93)],  where  <p  > is  the  average  value  of  p 

pGC  spec 


(p  ) = exp 

spec 


47th  sin  Y 


(3) 


One  might  conclude  from  this  formula  that  for  y sufficiently  small,  a 
surface  will  appear  smooth  ((p^)  - 1)  regardless  of  the  roughness. 
That  is,  the  Rayleigh  criterion  for  a smooth  surface  holds  for  arbitrary 
h.  We  shall  see  below  t!  at  since  shadowing  effects  have  thus  far  been 
neglected,  Eq.  (3)  will  not  generally  hold  for  any  h > 0 because  shadow- 
ing of  one  part  of  a random  rough  surface  by  another  tends  to  reduce  p 


^ • Quasispecular  Scattering  Theory 


This  statistical  model  postulates  a gently  undulating  random 
surface  where  the  radius  of  curvature  is  everywhere  much  greater  than  X. 
The  analysis  assumes  that  all  of  the  scattering  from  such 'a  rough  surface 
comes  from  "specular  points"  that  are  locally  align,,  such  that  they  re- 
flect signals  specularly  to  the  receiver.  This  model,  sometimes  known 
ns  the  cracked-egg  model  (boiled  egg),  is  analogous  to  moonlight  reflect- 
ing off  a slightly  rough  sea.  For  a Gaussian  distribution  of  surface 
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heights  and  a Gaussian  autocorrelation  function,  the  bistatic  radar  cross 
section  becomes  (Tyler  and  Ingalls,  1971) 


4 2 

q cot  p 


CT  (9  , 9,0)  = 

1 s s 


exp 


2 q 


2_ 
tan  9 

2 

2 tan  P 


R(y  , e ) 


(4) 


where  q = (k  - k )/k;  k and  k are  the  incident  and  scattered  wave 
is  1 s 

2 

vectors  in  Figure  3 and  k = 2rt/A.;  tan  |3  is  the  mean  square  unidirectional 
2 

slope  = (h/£)  ; h is  the  rms  height;  & is  the  correlation  length;  and 
2 2 2 

tan  P = q /q  , where  q and  q are  the  components  of  q parallel  and 

II  z II  z 

perpendicular  to  the  mean  surface  at  the  point  of  incidence.  |R(Y,e)|  is 
the  Fresnel  reflection  coefficient  for  specular  reflection  from  a perfectly 
smooth  dielectric  surface  (see  Barrick,  1970,  p.  700).  For  a perfectly 
conducting,  smooth  surface,  | R(y, e ) | -•  1,  The  scattered  signal  predicted 
by  this  model  has  both  coherent  and  noncoherent  components,  with  the  non- 
coherent component  dominating  as  g increases  (see  Beckmann  and  Spizzichino, 
1963,  Chapter  5 and  Figure  4 above). 

Quasi  specular  theory  has  been  developed  over  a number  of  years, 
bv  Beckmann  and  Spizzichino  (1963),  who  refer  to  quasispecular  as  diffuse 
scattering,  and  by  others,  with  a number  of  variations — e.g.,  non-Gaussian 
statistics.  A large  variety  of  lunar  and  terrestrial  radar  data  are  in 
reasonably  good  agreement  with  quasi  c-ulai  theory  as  shown  by  Barrick 
(1970,  pp.  753).  However,  quasispecu  ieory  breaks  down  at  low  grazing 

angles  unless  shadowing,  and  possibly  other  effects,  are  included.  A 
comparison  of  quasispecular  theory  including  shadowing  with  relevant  lunar 
bistatic  data  is  made  later  in  this  section.  For  certain  geometries  and 
surfaces  a "diffuse  component"  due  to  small-scale  roughness  may  dominate 
the  quasispecular  component,  as  discussed  below. 
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4. 


Shadowing  Theory 


At  small  grazing  angles  the  shadowing  of  one  part,  of  a surface 
by  another  and  multiple  reflections  between  surface  elements  become  in- 
creasingly important  for  both  geometrical  and  statistical  type  scattering 
models.  In  the  case  of  statistical  models,  the  multiple-scattering  problem 
has  not  yet  been  attacked  quantitatively,  but  considerable  efforts  have 
been  made  on  the  shadowing  problem — e,g. , Beckmann  (1965),  Brockleman  and 
Hagfors  (1966),  Smith  (1967),  and  Bass  and  Fuchs  (1972,  in  Russian). 

The  chapters  on  shadowing  in  Bass  and  Fuchs  (1972)  have  recently 
been  translated  by  C.B.  Vesecky,  but  are  not  generally  available  in  English 
so  a brief  comment  on  them  is  appropriate.  They  begin  with  a general  for- 
mulation that  includes  both  shadowing  and  multiple  scattering,  but  special- 
izes to  single  scattering.  They  characterize  the  relevant  properties  of  a 
random  surface  with  normally  distributed  height  and  slope  by  a function  A 
as  follows: 


whe  re 


a 


tan 


tan 


Erf c (X) 


(6) 


and  tan  3 is  the  unidirectional  rms  surface  slope  for  a Gaussian  surface, 
o 

The  function  A is  also  defined  for  the  variable  b = tan  yg/  | /2  tan  3q) • 


IF 
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Bass  and  Fuchs  develop  a reflection  coefficient  (PSpec^  for  the  average 
reflected  field  in  the  limiting  case  of  weak  shadowing  (a  » 1)  and  near- 
specular  reflection 


<D  > = 

spec 


1 


Erfc  (vi ) 


[1  - 2A(a)]  e 


-2(kh  sin 


V 


(7) 


Note  that  Eq.  (7)  is  simply  Eq.  (3)  with  a multiplicative*  correction 
factor  to  account  for  shadowing.  In  the  limiting  case  of  strong  shadow- 
ing (a  « 1),  (Pspec**  becomes 


r 

2 

I exp 

- — + iq  hx-A  (a)Erfr (x) 

2 z 

(8) 


They  also  derive  results  for  the  average  powoi  reflection  coeffi- 
cient by  calculating  a function  Q(y  , y ) such  that  the  reflected  power 

i s 

with  shadowing  included  is  Just  Q times  the  power  reflected  without  shadow- 
ing. In  the  case  of  strong  shadowing  (a  « 1)  we  have  for  forward  scatter: 


Q(Y.,Y  ) 
i s 


1 - exp {-f A(a)  + A(b)] } 
A(a)  + A(b)  


(9) 
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and  for  backward  scatter: 


Q(Y . , Y ) 

1.  s 


-A  (c) 
1 - c 

A ( c ) 


(10) 


where  c is  the  minimum  of  a and  b. 

Following  these  results,  Bass  and  Fuchs  examine  the  case  where 
a « 1,  and  they  derive  effective  distribution  functions  for  the  height 
and  slope  based  on  the  work  of  Smith  (1967).  A new  shadowing  function  Q 
is  derived  for  forward  scatter: 

^(Yi’  V = 1 + A(a)  + A (b)  ( 

and  for  monostatic  backscatter: 


Q(Yt>  = 


1 

1 + A (a) 


(12) 


The  backscatter  result,  Eq.  (12),  is  in  agreement  with  Smith  (1967).  A 
comparison  of  the  effective  probability  distributions  associated  with 
Q and  Q and  the  simulation  experiment  of  Brockleman  and  Hagfors  (1966) 
favors  the  choice  of  §.  However,  the  choice  does  not  seem  to  be  critical 
for  the  cases  compared.  Clearly,  this  Soviet  effort  on  the  shadowing  prob- 
lem is  substantial,  and  the  authors  arc  fully  aware  of  published  Western 
results.  However,  the  final  results  are  essentially  the  same  as  those  of 
Smith  (1967)  in  the  case  of  intermediate  (a  « 1)  and  strong  (a  « 1) 
shadowing.  The  fact  that  they  do  not  compare  their  theory  with  Soviet 
experiments  probably  indicates  that  such  experiments  do  not  exist  as  far 
as  the  open  Soviet  literature  is  concerned. 
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Needless  to  say,  shadowing  is  extremely  important  for  low-angle 
radar  tracking.  Later  in  this  section  we  compare  this  shadowing  theory 
with  forward-scatter  data  collected  during  an  Apollo  16  bistatic  radar 
experiment  that  explored  the  lunar  surface.  We  will  see  that  for  grazing 
angles  ^ r.nd  V g)  smaller  than  or  comparable  to  the  rms  surface  slope 
(i.e.,  a < 1)  scattering  by  the  quasispecular  process  is  significantly 
reduced  by  shadowing.  This  is  especially  important  because  one  would 
expect  an  increase  in  the  bistatic  radar  cross  section  as  y and  y ap- 
proach zero,  on  the  basis  of  quasispecular  theory  for  a plane  surface 
without  shadowing.  For  a sufficiently  curved  surface,  the  divergence 
factor,  D,  reduces  the  scattered  power  seen  by  an  observer  at  low  grazing 
angles — see  Section  Il-C  below. 

5.  Diffuse-Scattering  Theory 

The  quasispecular  scattering  from  a rough  surface  described 

above  applies  to  scattering  angles  near  the  specular  direction  (9  0 ) 

s ~ i 

In  addition,  there  is  usually  a diffuse  scattering  process  due  to  rough- 
ness on  a scale  \,  a scale  much  smaller  than  the  roughness  giving  rise 
to  the  quasispecular  component.  This  diffuse  component  will  generally 
dominate  the  quasispecular  component  for  scattering  angles  far  enough  away 
from  the  specular  direction.  At  present  there  is  no  widely  accepted  sta- 
tistical type  analysis  of  this  diffuse  component  and  one  is  forced  to  fall 
back  on  a semi-empirical  model  such  as  the  generalized  Lambertian  law 


o „ n, 

CT  = K (cos  9 cos  9 ) 

1 i s 


(13) 
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or  the  Lommel-Seeliger  law 


CT°  = K 


(2  cos  9 cos  9 \ 

i £) 

,■  )s  9 + cos  9 / 

^ i s / 


(14) 


K , n and  are  free  parameters  adjusted  to  fit  a given  experimental 
data  set.  Bar.-  tck  (1970,  p.  677)  comments  that  these  two  laws  compare 
favorably  with  some  backscatter  data  in  directions  where  diffuse  scatter- 
ing is  thought  to  dominate.  At  low  grazing  angles,  Barton  (1974,  p.  690) 

o 

uses  what  is  equivalent  to  the  Lambertian  law  with  K = 4jt(h/P  X)  and 

1 o 

n,  = 1 in  his  analysis  of  diffuse  scattering  in  the  low-angle  radar  track- 
ing problem. 

6.  Scattering  from  Vegetation 

As  one  can  well  imagine,  vegetation  is  difficult  to  model  as 
a radio-wave  scatterer.  A successful  model  must  account  for  several 
general  characteristics  observed  experimentally,  as  follows.  The  inci- 
dent wave  is  attenuated  as  it  propagates  into  the  vegetation  layer  so 
that  the  depth  of  the  layer  is  not  important  so  long  as  it  is  greater 
than  \ (Peake,  3959a).  (This  implies  that  a vegetation  layer  does  not 
have  to  be  very  thick  to  effectively  mask  the  underlying  terrain.)  The 
scattering  properties  must  also  depend  on  the  dielectric  properties, 
moisture  content,  size,  and  density  of  the  vegetation.  Since  vegetation 
is  often  of  a linear  form  (e.g.  , grass,  pine  needles,  etc.),  one  expects 
some  significant  dependence  of  the  return  on  the  polarization  of  the  in- 
cident wave. 

Barrick  (1970,  p.  689)  reviews  a model  by  Peake  (1959a  and  1959b) 
in  which  vegetation  is  modeled  by  a collection  of  thin  dielectric  rods, 
randomly  arranged,  but  preferring  the  vertical,  and  terminated  on  the 


uppor  end  by  a flat  horizontal  plane.  The  model  falls  in  the  geometric 
class  of  models  and  is,  according  to  Barrick,  the  only  model  that  can 
qualitatively  and  quantitatively  describe  scattering  from  vegetation-covered 
surfaces.  As  with  the  quasispecular  theory  above,  it  does  not  take  account 
o^  shadowing  and  multiple  scattering  and  so  is  likely  to  fail  at  the  small 

glazing  angles  (V  in  Figure  3)  that  are  of  interest  in  the  low-angle  radar 
tracking  problem. 

7.  Components  of  the  Resultant  E Field  at  the  Radar  Receiver 

Let  us  think  of  the  electric-field  vector  at  the  receiver  R in 

Figure  2 as  the  sum  of  three  components:  The  E field  (E  ) backscattered 

from  the  target  along  TR,  the  E field  (E  ) coherently  reflected  at  the 

R 

specular  reflection  point  (where  0 = 9 v , and  the  scattered  E field  (E  ) 

S ' s 
which  is  incoherent.  At  a given  instani  m time  the  resultant  E would 

appear  somewhat  as  in  Figure  5,  depending  on  the  relative  magnitudes  of 
the  three  components.  We  will  assume  here  that  all  three  components  are 
received  with  approximately  the  same  antenna  gain  as  in  the  case  of  low- 
angle  radar  tracking.  If  the  surface  is  only  slightly  rough,  specular 
reflection  dominates,  Eg  ~ 0,  Er  ~ E^  and  the  resultant  E varies  between 
about  0 and  2 Et  according  to  the  path  difference  between  r and  rj  + r2 
in  1 igure  2.  Now  if,  in  this  case,  the  target  T moves,  E varies  systemati- 
cally from  ~0  to  ~2  Et  according  to  the  changes  in  path  difference,  with 
disastrous  results  if  a monopulse  tracking  radar  is  used  (see  Hey  and 

Parsons,  1955,  and  Evans,  1966).  For  the  rough  surface  case  E ~ 0 

R ’ 

ET>  V and  E varies  between  E + E and  E - E . In  this  case 
1 & i s T s 

— » — ► — * 

E = E^  + E^,  and  as  the  target  moves  (or  the  sea  surface  changes)  E 

s 

varies  in  a random  manner,  but  with  a certain  autocorrelation  time  T 

’ s 

i.e.,  a characteristic  time  during  which  the  scattering  surface  changes 
from  one  realization  of  an  ensemble  of  random  surfaces  to  another,  and 
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hence  during  which  E^  changes  significantly.  Beckmann  and  Spizzichino 
(1963)  consider  the  rough-surface  case  (p,  120)  and  estimate  (p„  282)  a 
horizontal  autocorrelation  distance  A*  x that  we  can  convert  to  an  auto- 
correlation time  t bv  dividing  by  the  target's  velocity  (V): 


where  h is  the  greater  of  the  two  heights  h^  and  h^  in  Figure  2.  The 
period  of  the  variation  T in  the  slightly  rough  (specular  reflection) 
case  (p.  283)  is  similarly 


Ax  2 

^ s r X 

T = sa 

R V 2 Vh  h 
r t 


(16) 


FIGURE  B COMPOSITION  OF  THE  ELECTRIC-FIELD  VECTOR  AS  SEEN  AT  THE 

RADAR,  R,  IN  FIG.  2.  The  electric  field  at  the  radar,  E,  is  the  sum  of 
— ► 

the  direct  wave,  ET,  from  the  target,  the  reflected  or  coherently  scattered 
wave,  Er,  and  the  incoherently  scattered  wave,  Es. 


Consider  the  case  of  a target  in  level  flight  over  land  where 

ht  = m’  hr  = ® m>  v = 250  m s (560  mph),  r = 20  km  and  X = 3 cm 

(X-band).  Thus,  from  Eqs,  (15)  and  (16)  we  have  T ~ 10  s and  T ~ 48  s 

s R 

If  we  accept  Eqs.  (15)  and  (16)  as  correct,  E^  varies  extremely  slowly 
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on  the  time  scale  of  a typical  radar  pulse  repetition  period  10  s). 

In  fact,  for  times  small  compared  to  T , specularly  reflected  (E  ) and 

s R 

scattered  (E  ) components  look  virtually  the  same  to  the  radar  receiver 
s 

except  for  magnitude.  Beckmann  and  Spizzichino  0963,  p.  284)  point  out 
that  T may  be  shorter  than  given  by  Eq.  (3)  due  to  atmospheric  scintilla- 
tion effects.  Indeed,  measurements  by  Bullington  (1954,  p.  1259)  indicate 

that  T should  be  about  0.5  s for  one  particular  case, 
s 

This  analysis  would  indicate  that  over  land,  care  must  be  taken 

in  averaging  signals  to  reduce  the  errors  caused  by  incoherent  scattering. 

This  is  because  the  scattered  field  at  the  receiver  will  not  vary  appreciably 

on  time  scales  shorter  than  T and  therefore  averaging  over  times  less  than 

s 

T will  not  be  effective  in  reducing  u-nwanted  random  signal  fluctuations. 

s 

8.  Doppler  Shift  of  the  Direct  and  Scattered  Echoes 

The  radar  echo  from  a rnovir-  target  is  Doppler-shifted  by  an 

amount  Af  = f - f where  f is  the  observed  frequency  and  f is  the  trans- 
° o 

mitted  frequency.  Referring  to  Figure  2 , Af  = (fQ/c)  (r  • V + d • V) 
where  c is  the  velocity  of  light,  r is  a unit  vector  along  r,  and  d is 
a unit  vector  along  the  direction  with  which  the  echo  signal  leaves  the 
target.  For  the  direct  echo  along  path  TR  we  have  5=5,  whereas  for  the 
scattered  echo  van  path  TOR,  d = (r  /|r  |).  Since  Af  will  vary  with  the 
direction  of  d,  it  is  in  principle  possible  to  distinguish  between  the 
direct  and  scattered  echoes  on  the  basis  of  Doppler  shift. 

To  appreciate  the  significance  of  Doppler  shift  in  the  low-angle 
tracking  problem,  we  will  consider  the  somewhat  simplified  case  illustrated 
in  Figure  6.  Here  we  have  assumed  a flat  earth  and  will  consider  only 
echoes  that  are  scattered  in  the  plane  of  the  figure.  In  spite  of  these 
limitations  some  interesting  features  will  be  evident.  Applying  the  for- 
mula for  Af  given  above  to  the  geometry  of  Figure  6 we  have,  for  the  direct 
echo  along  TR, 
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Af  = ( f V/c)(2  cos  o)  = (f  V/c) 
° o 


(17) 


and  for  the  scattered  echo  along  TOR, 


Af  = (f  V/c)(cos  a + cos  B) 
o r ' 


and 


Af  = (f  V/c) 

COS 

tan'1/  \ 

° 

L \ c l\ 

+ cos 

\G  - X / 

- \ s /J 

(18) 


a 


FIGURE  6 GEOMETRY  FOR  CALCULATING  THE  DOPPLER  SHIFT.  The  radar  echoes  arrive  , 
the  radar  R,  directly  from  the  target,  T.  via  path  TR  and  scattered  from  the  surface- 
e.g.,  via  the  path  TOR.  G is  the  ground  range  from  radar  to  target,  and  X is  the 
ground  range  to  a scattering  point,  0 The  mean  surface  of  the  earth  is  assumed  to 
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To  apply  these  formulae  to  a typical  low-angle  tracking  problem,  we  take 
a taiget  velocity  V = 250  ms  (560  mph) , a radar  wavelength  \ = 3 Cm,  a 
radar  height  = 5 m,  a target  height  ^ = 100  m,  and  target  ground  ranges 
G = 10  and  20  km.  In  Figure  7 we  have  plotted  the  Doppler  shift  of  the 
scattered  echoes  as  a function  of  Xg , the  ground  range  of  the  scattering 
point,  and  compared  it  to  the  Doppler  shift  of  the  direct  echo.  Consider- 
ing the  curvature  of  the  earth  and  standard  refraction  of  the  radar  signal, 
we  see  that  the  real  horizon  for  h^  = 5 m is  about  9 km.  If  we  consider 
scattering  points  above  the  mean  surface,  the  horizon  moves  to  larger 

ranges— e.g. , for  scattering  points  1 m above  the  mean  surface  the  real 
horizon  is  about  13  km. 

It  is  clear  from  Figure  7 that  Doppler  discrimination  between 
the  direct  and  scattered  echoes  becomes  increasingly  difficult  as  X de- 
creases. However,  Doppler  discrimination  could  well  be  useful  in  discrimi- 
nating against  echoes  scattered  near  the  target  (Xg  - G) . Echoes  scattered 
in  this  region,  known  as  the  "horizon  component"  of  the  scattered  echo, 
can  in  fact  be  relatively  important,  as  noted  in  connection  with  Eq.  (2). 
Such  discrimination  would  be  less  useful  for  large  values  of  G (e.g., 

20  km),  since  much  of  the  area  near  the  target  would  be  below  the  radar’s 
horizon  anyway.  Only  those  scattering  points  lying  7 m and  more  above  the 

mean  surf ace-e . g . , small  hills-would  be  visible  at  X « G = 20  km 

s * 

For  both  values  of  G considered  here  the  specular  reflection 
point  is  very  near  the  radar  (Xg  ^ 1 km),  so  for  a very  smooth  surface 
where  specular  reflection  is  important,  Doppler  discrimination  would  be 
relatively  difficult.  As  the  surface  becomes  rougher  the  echo  scattered 
quasispecularly  comes  from  a broader  range  of  values  of  Xg,  some  regions 
being  more  important  than  others-cf.  Barton  (1974,  p.  691) , and  Beckmann 
and  Spizzichino  (1963,  Chapter  12).  Since  a°  for  diffuse  scattering 
varies  more  slowly  with  position  than  the  factor  (1/r  r ) in  Fq.  (2),  the 
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echo  arising  irom  a diffuse  process  will  tend  to  be  concentrated  near  the 
horizon  and  near  the  radar.  In  general  the  scattered  echo  will  be  com- 
posed of  a sum  of  echoes  from  points  with  different  Doppler  shifts,  each 

bo i ng  weighted  according  to  the  size  of  the  bistatic  radar  cross  section 
o „ 

o . As  mentioned  above,  the  theoretical  models  for  o relevant  to  the 

low-angle  tracking  problem  are  not  at  all  well  established;  hence,  only 

qualitative  statements  are  probably  justified.  However,  once  adequate 

low-angle  scattering  models  become  available,  the  Doppler  profile  (echo 

intensity  vs.  frequency)  can  be  easily  calculated. 

Both  Barton  (1974,  p.  694)  and  Fjeldbo  (1964)  have  calculated 
Doppler  profiles  for  specific  scattering  laws.  Barton's  calculation  uses 
a scattering  law  that  approaches  the  Lambertian  law  [Eq.  (13)]  for  small 
grazing  angles  and  is  specifically  addressed  to  the  low-angle  tracking 
problem.  He  finds  the  half-power  bandwidth  Af  to  be  about  1 Hz  or  less, 
depending  on  the  specific  case.  Only  for  the  case  of  a very  rough  sur- 
face (hA)  » 20  and  a 105-m-al  ti tude  target  at  10  km  range  does  he  find 
appreciable  scattered  echo  power  more  than  2 Hz  from  the  direct  target 
echo.  In  considering  these  results  one  should  remember  that  the  corre- 
spondence between  reality  and  the  Lambertian  scattering  law  used  in  these 
calculations  is  not  well  known  at  low  angles.  However,  it  is  probably  as 
good  a choice  as  any  other  at  the  present  time. 

Fjeldbo  (1964)  derived  a formula  for  the  half-power  bandwidth  of 

the  scattered -echo  spectrum  using  a statistical  model  similar  to  the  quasi- 

specular  scattering  law  of  Eq.  (4).  He  assumes  V = y = y and  finds 

i s 

Af  = 4 y^2  In  2 (V  /A.)  tan  (8  ) sin  y (19) 

so 

where  is  the  velocity  of  the  specular  point  on  the  mean  surface.  Using 

the  geometry  of  Figure  6 we  find  that  V = [h  /(h  + h )]  V where  V is 

s r r t 
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assumed  parallel  to  the  surface.  Using  the  same  parameters  as  for  Figure  7 
(i.e.,  X = 3 cm,  G 10  km,  = 5 m,  ht  = 100  m,  and  V = 250  ms"1),  and 
taking  tan  ^ = 0.04  (corresponding  to  moderately  rough  terrain),  we  have 

Yi  = 0,01  rad  ar‘d  ^ = °'75  Hz*  HorG  again  we  do  not  know  how  realistic 
Fjeldbo's  scattering  law  is  at  low  angles.  Fjeldbo's  result  has  been  used 
in  connection  with  lunar  bistatic  radar  data  (Tyler  and  Simpson,  1970)  and 
predicts  reasonable  values  of  the  unidirectional  rms  slope  (tan  3 ) for 
Y.  and  Yg  » 8q.  However,  at  small  grazing  angles  the  theory  appears  to 

break  down,  probably  because  shadowing  was  neglected  (see  Tyler  and  Ingalls, 
1971,  p.  4775). 

Cased  on  these  analyses  one  would  have  to  use  very  narrowband 
filters  («?  1 Hz)  to  discriminate  between  direct  and  scattered  echoes  by 
means  of  Doppler  shift.  But  if,  as  one  might  expect  from  quasispecular 
scattering  theory,  the  scattered  echo  for  low-angle  targets  turns  out  to 
come  from  a relatively  small  number  of  unshadowed  "specular  points,"  each 
one  would  have  a more  or  less  unique  Doppler  shift  that  would  change  as 
the  target  moves.  It  seems  possible  that  one  could  exploit  this  situation 
by  placing  the  radar  at  a site  such  that  the  Doppler  shifts  of  the  domi- 
nant "specular  points"  were  as  far  as  possible  from  the  Doppler  shift  of 
the  target.  Alternatively,  one  might  observe  the  "Doppler  tracks"  of 
known  specular  points"  or  even  of  transponders  to  obtain  information  con- 
cerning target  location.  Some  really  definitive  experimental  work  on 
scattering  at  small  grazing  angles  is  necessary  to  evaluate  such  a possi- 
bility. Some  suggestions  regarding  experimental  work  are  given  in  Section 
IV  below. 

9*  Deficiencies  in  Scattering  Theory  Relevant  to  Low-Angle 

Radar  Tracking 

The  outstanding  problems  in  scattering  from  rough  surfaces  at 
low  grazing  angles  are  as  follows: 
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(1)  Qua  si  specular  scattering  theory  breaks  down  at  low 
grazing  angles  because  shadowing  and  multiple  scattering 
have  been  neglected.  Some  theoretical  work  on  shadowing 
exists,  but  is  untested  by  experiment  (except  as  done  for 
the  lunar  surface  in  Section  II-C-4);  and  multip:e  scatter- 
ing is  still  neglected. 

(2)  No  good  statistical  or  geometric  theory  exists  for  the 
diffuse-scattering  component  for  any  grazing  angle.  To 
the  authors'  knowledge,  semi-empirical  models  such  as 
Eqs.  (13)  and  (14)  have  not  been  compared  with  experi- 
mental data  at  low  grazing  angles. 

(3)  Scattering  from  vegetation  is  clearly  a difficult  problem, 
especially  at  low  grazing  angles.  Peake's  model  (Barrick, 
1970,  p.  689)  appears  to  be  the  only  model  known  that  can 
"describe  qualitatively  and  quantitatively  the  scattering 
by  vegetation-covered  surfaces."  But  Peake's  model  de- 
viates from  measurements  at  low  grazing  angles. 

Once  scattering  laws  for  a particular  range  of  parameters  are 
well  understood  from  both  theoretical  and  experimental  viewpoints,  one 
can  predict  with  confidence  how  the  scattered  echo  will  differ  from  the 
direct  echo--e.g.,  in  time  delay,  Doppler  shift,  intensity,  polarization, 
etc. — and  algorithms  to  do  the  discrimination  can  be  derived  and  evaluated. 
At  present  the  above  deficiencies  in  scattering  theory  will  be  reflected 
as  uncertainties  in  predicted  tracking  errors  for  a given  radar  system. 


C*  Expei-imotita1  Data  on  Microwave  Scattering  from  Reugh  Surfaces. 
and  Comparison  with  Theory  ~ " ' ’ 

1.  Ge  ne  ra  1 

Experiments  gathering  data  on  the  scattering  of  microwave  radia- 
tion irom  the  Earth's  land  and  sea  surfaces  have  been  going  on  for  consider- 
ably more  than  20  years.  In  addition,  radar  waves  have  been  used  to  explore 
+he  Moon's  surface  and  the  surfaces  of  Mars  and  Venus.  This  mass  of  ex- 
perimental data  has  encouraged  theorists  to  investigate  the  problem  of 
electromagnetic-wave  scattering  from  randomly  rough  surfaces.  There  have 
been  notable  cases  of  a theory  successfully  explaining  a set-  of  experi- 
mental data— e.g.,  quasi  specular  scattering  theory  and  lunar  radar  data. 
However,  when  one  attempts  to  compare  theory  with  available  experimental 
data,  there  are  a number  of  inherent  problems. 

First,  experimental  results  are  often  given  in  terms  of  the 
power  reflection  coefficient  (p2)  observed  under  a particular  set  of 
experimental  conditions.  As  shown  in  Section  II-C-4  below  and  in  Barton 
(1974),  a theoretical  calculation  of  p2  involves  an  integration  over  the 
scattering  surface  as  well  as  specialization  of  the  theory  to  match  the 
experimental  conditions.  So  only  one  implication  of  the  theory  is  tested 
and  the  point  of  comparison  is  an  integral  quantity.  Such  a comparison 
is  clearly  not  an  ideal  one;  yet  it  is  often  all  that  is  possible. 

Another  difficulty  arises  in  the  experimental  measurement  of 
the  scattered  wave.  A full  specification  of  the  scattered  wave  requires 
the  measurement  of  the  coherency  matrix  [J]  (or  the  general  Stokes  param- 
eters) as  a function  of  direction  of  arrival,  time  delay,  and  Doppler  shift 
(see  Born  and  Wolf,  1975,  pp.  544-555,  and  Hagfors,  1967).  Experimental 
measurement  of  such  a full  specification  is  a substantial  task  and  is 
generally  attempted  only  in  radar  astronomy  experiments  such  as  that  by 
Tyler  and  Howard  (1973)  discussed  below.  Thus  the  experimental  data  that 
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would  enable  one  to  make  a crucial  comparison  of  theory  with  experiment 
are  often  not  available  at  all  or  only  with  substantial  data  processing. 

A genuine  test  of  theory  by  experiment  must  include  a sufficiently 
accurate  characterization  of  the  scattering  surface.  For  example,  the  quasi- 
specular  scattering  theory  developed  by  Beckmann  in  Beckmann  and  Spizzichino 
(1963)  and  others  requires  the  unidirectional  rms  slope  of  the  scattering 
surface  (tan  P ) on  a scale  large  compared  with  the  wavelength — say,  about 
10  A.  . A natural  surface  will  also  possess  roughness  on  scales  shorter 
than  X as  well  as  on  scales  very  much  longer  than  X — e.g.,  hills,  etc. 

At  wavelengths  in  the  VHF  and  UHF  ranges  even  subsurface  characteristics 
can  be  important.  So  in  a field  experiment  one  must  be  careful  to  measure 
the  right  parameters  over  the  right  scale  lengths.  Only  then  can  meaning- 
ful comparison  between  theory  and  experiment  be  made. 

An  incisive  experiment  that  will  genuinely  test  a scattering 
theory  must  involve  a very  careful  choice  of  the  parameters  to  be  measured, 
both  of  the  scattered  radiation  and  of  the  scattering  surface.  While  it 
would  be  expensive  and  probably  unnecessary  to  measure  all  the  conceivable 
parameters  of  interest,  it  should  not  be  beyond  the  wit  of  man  to  design 
clever  experiments  that  establish  the  range  of  parameters  over  which  the 
critical  assumptions  of  a theory  are  valid.  In  fact,  a good  experiment 
that  establishes  the  weak  points  of  a particular  theory  could  very  well 
point  to  the  faulty  assumptions  that  need  to  be  rectified  in  subsequent 
theoretical  developments.  In  the  case  of  scattering  from  rough  surfaces 
at  very  low  grazing  angles,  the  critical  experiments  needed  to  test  current 
theories  have  yet  to  be  made,  so  far  as  we  have  been  able  to  determine. 

Some  specific  experimental  recommendations  are  made  in  Section  IV. 

2 . Data  on  Forward  Scattering  over  Land 

It  is  beyond  the  scope  of  this  report  to  survey  the  literature 
on  the  forward  scattering  of  microwaves  over  land.  However,  it  is  useful 
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to  cite  some  examples  of  relevant  work  and  to  comment  on  their  strong 
points  ar.ci  weaknesses.  One  set  of  experimental  data  used  by  Barton  (1974) 
for  comparison  with  theory  is  that  of  Bullington  (1954).  This  data  set  is 
really  a byproduct  of  the  route  survey  for  the  transcontinental  microwave 
relay  system.  Certainly  the  data  are  useful,  especially  since  a large 
number  of  paths  were  measured.  However,  these  data  cannot  be  used  in  a 
critical  test  of  theory  because  too  many  things  apparently  were  not  mea- 
sured. For  example,  the  nature  of  the  scattering  surface  is  recorded  only 
on  horizontal  scales  longer  than  about  1/2  mile.  Bullington  in  fact  con- 
cludes that  "...the  magnitude  of  the  reflection  coefficient  cannot  be 
predicted  accurately  from  the  gross  features  of  the  path  profile."  Of 
course,  these  experiments  were  never  designed  to  test  scattering  theory 
and  one  should  be  grateful  that  they  were  published  at  all.  Even  so,  they 
are  oiten  quoted  in  comparisons  of  theory  and  experiment. 

Barton  (1974)  quotes  reflection-coefficient  measurements  made 
by  McGavin  and  Maloney  (1959)  over  dry  Colorado  range  land  at  about  1 GHz. 
He  then  carries  out  a theoretical  calculation  that  yields  power  reflection 
coefficients  of  0.16  to  0.25  and  compares  these  to  the  mean  value  of  0.12 
for  the  experimental  data.  This,  Barton  (1974)  concludes,  "constitutes 
good  agreement,  considering  the  uncertainties  in  antenna  patterns,  surface 
slopes,  masking  and  illuminating  regions."  Certainly  the  theoretical  and 
experimental  values  are  reassuringly  close.  However,  such  a comparison 
oi  reflection  coefficients,  which  are  integral  quantities,  is  far  from  a 
crucial  test  of  the  theory. 

**•  Data  on  Forward  Scattering  over  the  Sea 

Here  again  our  purpose  is  only  to  comment  on  a few  examples  of 
this  typo  oi  data,  rather  than  attempt  any  survey.  The  experiment  sum- 
marized by  Beard  (1961)  shows  many  good  features  in  terms  of  comparing 
theory  with  experimental  measurement.  A number  of  properties  of  the 
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scattered  radiation  were  measured,  including  a division  into  coherent 
and  incoherent  portions.  The  scattering  surface  was  the  ocean  and 
various  wave  parameters  were  met  .'red  so  as  to  characterize  the  surface 
by  the  "apparent  ocean  roughness  parameter  (hyA ) , where  h = rms  wave 
height,  y = grazing  angle,  and  \ = wavelength.  Thus  the  experimental  data, 
gathered  on  transmissions  between  oil  drilling  platforms  in  the  Gulf  of 
Mexico,  could  readily  be  compared  +o  the  theoretical  predictions  available 
in  1961.  It  is  worthwhile  to  point  out  that  the  experiments  were  planned 
by  one  laboratory  but  executed  by  another  organization  in  close  coopera- 
tion. Possible  future  experiments  might  well  benefit  from  the  combined 
talents  of  two  organizations. 

The  measurement  of  both  coherent  and  incoherent  components  of 
the  scattered  wave  was  an  attempt  to  isolate  the  radiation  scattered  by 
a diffuse  scattering  mechanism.  Clearly  this  was  a good  idea  for  ap- 
plication to  low-angle  radar  tracking. 

From  the  viewpoint  of  more  modern  ideas  in  scattering  theory, 
these  Gulf  of  Mexico  experiments  were  inadequate  in  a number  of  respects. 
However,  our  comments  here  should  not  be  taken  as  criticism  of  the  original 
work;  but  rather  as  ideas  as  to  how  a future  experiment  along  similar  lines 
might  be  made  more  pertinent  to  current  scattering  theory.  Concerning  the 
scattered-wave  measurements,  a more  nearly  complete  characterization  of 
the  scattered  wave  is  needed.  For  example,  Tyler  and  Howard  (1973)  mea- 
sured the  deterministically  polarized  and  randomly  polarized  portions  of 
electromagnetic  waves  forward-scattered  from  the  lunar  surface.  This 
measurement  was  an  attempt  to  separate  radiation  scattered  by  a "quasi- 
specular  process,  which  is  relatively  well  understood,  from  that  scattered 
by  a diffuse  process,  which  is  not  well  understood  (see  Section  II-B). 
Also,  the  characteristics  of  the  ocean  surface  need  to  be  moir  carefully 
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measured.  For  example,  the  rms  wav-  height  was  measured  in  the  Gulf  of 
Mexico  experiment,  while  the  rms  wave  slope  is  the  quantity  needed  to 
relate  to  quasi  specula r theory.  The  characterization  of  the  surface 
roughness  also  needs  to  be  specified  over  a range  of  horizontal  distance 
scales  both  larger  and  smaller  than  the  wavelength. 

Lunar  Histatic  Radar  Experiment 

Tyler  and  Howard  (1973)  investigated  the  lunai  surface  using 
s-boml  and  VHF  transmitters  on  the  Apollo  (11,  15,  and  16)  command  mod- 
ules and  Earth-based  receivers  In  a blstatic,  CW  radar  experiment.  Htc 
experiment  geometry  is  shown  in  Figure  8.  They  have  kindly  provided 
some  of  the  Apollo  16  s-ba„d  data  for  use  in  this  report.  One  of  the 
unique  features  of  the  experiment  was  the  measurement  of  the  coherency 
matrix  of  the  scattered  radiation.  This  was  accomplished  In  practice 
by  recording  the  complex  components  of  both  the  right-  and  the  loft- 
circularly  polarised  portions  of  the  scattered  radiation.  From  this 
basic  measurement  the  scattered  signal  could  bo  resolved  Into  a deter- 
ministically polarized  (or  "polarized")  portion  associated  with  the 
quasispccular  scattering  mechanism  and  a randomly  polarized  (or  "un- 
polarized  ) portion  associated  with  ono  or  more  "diffuse"  scattering 
mechanisms.  When  adjusted  to  normal  incidence,  the  ratio  of  "polarized" 
to  "unpolarized"  power  was  between  5 and  9 dB  at  S-band.  Thus,  quasi- 
specular  scattering  is  the  dominant  mechanism  at  this  frequency.  inter- 
preting the  "polarized"  component  of  scattered  radiation  as  being  due  to 
the  quasispccular  scattering  mechanism,  the  experimenters  were  able  to 
determine  the  rms  slope  and  dielectric  constant  of  the  near-specular 

region  as  it  moved  across  the  lunar  surface  with  the  motion  of  the  com- 
mand  modulo. 

Lunar  slope  measurements  based  on  bistatic  radar  observations 
at  2.2  „ wavelength  have  been  compared  with  photographic  results  (Tyler 
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SPECULAR  REFLECTION  POINT 


FIGURE  8 SCATTERING  GEOMETRY  FOR  APOLLO  16  BISTATIC  RADAR  EXPERIMENT 
For  the  moon,  R = 1738  km. 


et  al.,  1971).  The  photographic  and  bistatic  radar  methods  are  in  agree 
ment  provided  the  bistatic  radar  scale  length  L can  be  adjusted  over  a 
range  of  about  + 20%.  This  successful  "ground  truth"  verification  of 
bistatic  radar  methods  lends  credibility  to  the  analysis  which  follows. 


Using  the  polarized"  component  of  the  scattered  signal  and 
knowing  the  experimental  conditions,  lyier  and  Howard  (1973)  inferred 
a power  reflectivity  (p  ) for  the  lunar  surface  by  assuming  a quasispecu- 
lar  scattering  mechanism.  This  inferred  reflectivity  is  the  quantity 
used  ^comparison  with  theory  later  in  this  section  (Figure,  9 and  10). 
gh  p is  known  only  to  within  a multiplicative  constant,  we  will  not 

find  this  a serious  drawback  because  the  data  cover  a wide  range  of  graz- 
ing angles. 


This  bistatic  lunar  radar  data  is  of  particular  interest  here 
because  it  can  be  compared  to  the  shadowing  theory  of  Section  TI-B-4  with 
a relative  minimum  of  uncertainty.  In  addition,  the  data  cover  a range 
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of  grazing  angles  from  80° 
lunai  data  with  scattering 
Rather,  our  purpose  here  is 
at  small  grazing  angles  and 
work, 


to  2.5  . Clearly,  a detailed  comparison  of 
theory  is  beyond  the  scope  of  this  report, 
to  show  the  importance  of  shadowing  effects 
to  suggest  approaches  for  more  definitive 


The  theoretical  value  of  the  power  reflection  coefficient  is 
taken  to  be  of  the  form 
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where  C is  a normalizing  constant,  D2  is  the  divergence  factor,  Q is  the 
shadowing  function,  and  £(R2  + R2)  is  the  classical  Fresnel  reflection 
coefficient  for  power.  For  the  geometries  of  interest  here,  Tyler  and 
Ingalls  (1971)  have  shown  that  to  first  order  a quasi  specular  scattering 
law  gives  results  equivalent  to  the  Fresnel  coefficient  used  in  Fq.  (20) 
above.  Hence  we  have  used  the  simpler  Fresnel  expression.  The  methods 
by  which  each  of  the  factors  was  calculated  are  described  briefly  below. 


When  a wave  scatters  from  a convex  spherical  surface  as  opposed 
to  a plane  surface,  the  curvature  of  the  surface  causes  the  paths  of  re- 
flected rays  to  diverge.  Thus,  an  observer  receives  less  scattered  power 
from  a convex  surface  than  he  would  from  an  equivalent  plane  surface. 

This  divergence  effect  is  taken  into  account  by  the  divergence  factor  D2 
in  Eq.  (20).  The  geometry  is  illustrated  in  Figure  8,  where  we  have  as- 
sumed specular  scattering.  Kerr  et  al.  (1951,  p.  406)  quote  a result  of 
Van  do  Pol  and  Bremmer  (19.39)  giving 
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Given  the  experimental  geometry,  application  of  the  law  of  sines  yields 

2 

all  the  parameters  to  find  D from  Eq . (21). 

2 

The  inferred  reflectivity  p measured  by  Tyler  and  Howard  (1973) 

has  been  adjusted  to  take  the  Moon's  spherical  shape  into  account  and  thus 
2 

their  values  of  p given  in  Figures  9 and  10  refer  to  a scattering  surface 

which  is,  on  the  average,  flat.  Henc^  , in  calculating  theoretical  values 
2 

of  p to  compare  with  the  bistatic  radar  measurements,  we  have  used  the 

2 

divergence  factor  appropriate  to  a flat  surface  (D  = 1)  in  Eq.  (20). 


The  shadowing  function  Q takes  account  of  the  fact  that  one 

portion  of  a rough  surface  can  shadow  another,  especially  at  low  grazing 

angles.  A knowledge  of  the  grazing  angle  (y  = V - v ) and  the  unidirec- 

i s 

tional  rms  slope  (tan  9 ) determines  Q via  Eqs,  (6)  and  (11). 

o 


The  classical  Fresnel  reflection  coefficients  for  electric-field 
components  perpendicular  to  and  parallel  to  the  plane  of  incidence  are, 
respectively, 


_ sin  v - (e  - cos  ')' 

RJ.  2 "■ 

sin  y + (e  - cos  y)‘ 


(22) 


and 


R„  = 


e sin  v - (e  - cos  Y)‘ 


e sin  Y + (e  - cos  y)‘ 


(23) 


where  e is  the  average  dielectric  constant.  To  obtain  the  classical 
reflection  coefficient  for  power,  we  need  to  know  the  polarization  of 
the  incident  wave.  In  the  case  examined  here  the  incident  wave  is  ap- 
proximately equally  divided  between  parallel  and  perpendicular  components, 

i 2 2 

so  we  use  £(R^  + R ).  The  average  dielectric  constant  e is  assumed  to 
be  3 in  all  cases  for  the  lunar  surface.  Values  between  2.8  and  3.1  have 
been  observed,  however. 
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INFERRED  REFLECTIVITY,  p* — arbitrary  units 


FIGURE  9 COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  VALUES  FOR  THE 
POWER  REFLECTIVITY  (p2)  OF  THE  LUNAR  SURFACE.  The  ordinate  scale 
is  arbitrary  since  the  experimental  value  of  p2  is  known  only  to  within  a 
multiplicative  constant.  Tan  (30  is  the  unidirectional  rms  slope  of  the  lunar 
surface  assumed  for  the  theoretical  calculations  of  p2.  The  Apollo  16  data 

( line)  were  kindly  provided  by  Dr.  G.  L.  Tyler  of  the  Stanford  Center 

for  Radar  Astronomy. 


POWER  REFLECTIVITY.  pZ  — arbitrary  unit! 


GRAZING  ANGLE,  y — deg 


FIGURE  10  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  VALUES  FOR  THE 
POWER  REFLECTIVITY  <p2,  0F  THE  LUNAR  SURFACE.  The  ordinate  scale 
is  arbitrary,  since  the  experimental  value  of  p 2 is  known  only  to  within  a 
multiplicative  constant.  Tan  0O  is  the  unidirectional  rms  slope  of  the  lunar 
surface  assumed  for  the  theoretical  calculations  of  p?.  This  semi-log  presentation 
for  grazing  angles  between  0°  and  7°  allows  a closer  examination  of  the  role  of 

shadowing  at  small  grazing  angles.  The  Apollo  16  S-band  data  ( line)  were 

kindly  provided  by  Dr.  G.  L.  Tyler  of  the  Stinford  Center  for  Radar  Astronomy. 
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The  theoretical  values  of  p ' , bn  ,ed  on  Eq.  (20),  and  the  Apollo 

16  bistatic  radar  measurements  of  p are  compared  in  Figures  y and  10. 

Since  the  experimental  values  of  p2  are  known  only  to  within  a mulitpli- 

cative  constant  we  cannot  compare  the  theoretical  and  experimental  values 

on  an  absolute  basis.  However,  by  setting  the  normalizing  constant  C 

such  that  the  theoretical  curves  match  the  experimental  curve  at  near 

normal  incidence  (V  - 80°),  the  shapes  of  the  two  curves  can  be  compared. 

This  particular  normalization  is  chosen  so  as  to  make  the  experimental 

and  theoretical  curves  agree  at  large  grazing  angles  where  shadowing  ef-^ 

fects  are  unimportant.  No  significance  should  be  attached  to  the  value 

of  C since  the  vertical  scale  in  Figures  9 and  10  is  an  arbitrary  one. 

2 

To  use  Eq.  (20)  elsewhere,  as  a model  for  the  calculation  of  P in  an 
absolute  sense,  C should  be  taken  as  unity. 


The  three  theoretical  curves  differ  only  in  the  value  of  rms 

unidirectional  slope  (tan  0°)  used  to  calculate  the  shadowing  factor  Q. 

The  curve  for  0 =0°  includes  no  shadowing  of  fects  — i.e. , Q = 1 in 

o 

Eq.  (20).  It  allows  one  to  judge  the  relative  importance  of  the  shadow- 
ing function  Q.  On  the  Moon  values  of  0 range  fjom  about  2 to  4 for 
the  rolativei-..  smooth  maria  to  about  4°  to  8-  for  the  relatively  rough 
highland  regions.  Thus  the  curves  for  0q  = 4°  and  8 plotted  in  Figures 
9 and  10  correspond  to  rather  rough  mare  and  highland  terrain  respectively. 

For  comparison  the  value  of  0q  fpr  the  Sierra  Nevada  mountains  in  California 

o 

is  roughly  7 . 


2 

The  scattering  region  to  which  the  measurements  of  P apply 

moves  across  tlx?  lunar  surface  with  the  movement  of  the  transmitter  in 

the  command  module.  Each  point  along  the  track  of  the  scattering  region 

corresponds  to  a different  value  of  y because  the  experimental  geometry 

2 

changes  with  transmitter  movement.  Thus  measurements  of  P at  different 
values  of  Y in  Figures  9 and  10  correspond  to  different  scattering  regions 
on  the  lunar  surf  ace.' ~ For  high  values  of  y the  scattering  region  is 
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located  near  the  center  of  the  lunar  disk  passing  near  the  craters  Descartes 
and  Ptolcmaeus  as  y moves  to  lower  values.  As  y continues  to  decrease,  the 
scattering  region  moves  across  the  southern  end  of  Oceanus  Procellarum  and 
finally  over  the  western  limb.  For  low  grazing  angles,  which  are  of  prin- 
cipal interest  here,  the  scattering  region  is  the  lunar  highlands  near  the 
western  limb.  Hence  a value  of  g from  6 to  8 is  appropriate  to  the  in- 
terpretation of  Figures  9 and  10  for  y £ io° . 

At  large  grazing  angles  (y  2*  15°)  agreement  between  the  experi- 
mental and  theoretical  curves  of  Figure  9 is  excellent.  Below  y = 15° 
the  experimental  curve  for  g^  = 8°  agrees  most  closely  with  the  experi- 
mental data  and  the  theoretical  curve  without  shadowing  (g  =0°)  becomes 

o 

obviously  incorrect.  The  fact  that  the  g = 8°  curve  is  a better  fit  than 

o 

the  curve  for  g^  =4  is  expected  since  for  y 10°  the  scattering  region 
on  the  lunar  surface  is  a relatively  rough  highland  region  where  one  ex- 
pects to  find  g^  in  the  6°  to  8°  range. 

At  small  grazing  angles  (y  < 10°)  the  experimental  curve  falls 
progressively  further  below  the  theoretical  curves.  In  Figure  10  this 
feature  is  shown  in  more  detail.  The  main  difference  between  the  theore- 
tical curve  for  g^  = 8°  and  the  experimental  curve  in  Figure  9 is  that 
the  experimental  curve  breaks  downward  more  sharply  and  at  a larger  value 
of  y than  the  theoretical  one.  The  most  probable  explanation  for  this 
phenomenon  appears  to  be  the  neglect  of  a particular  type  of  multiple 
scattering  in  the  theoretical  model  of  Eq.  (20).  In  the  model  of  Eq.  (20) 
waves  that  are  scattered  toward  the  receiver  by  a properly  tilted  facet 
arc  only  shadowed  to  the  extent  thct  a random  fluctuation  from  a mean  plane 
surface  may  intercept  them.  That  is,  the  shadowing  model  docs  not  include 
the  effect  of  surface  curvature.  Thus  waves  that  scatter  from  properly 
tilted  facets  on  the  transmitter  side  of  the  specular  point  and  that 
escape  shadowing  by  random  surface  fluctuations  are  assumed  to  propagate 
to  the  receiver.  In  reality  the  raj  path  of  such  a wave  may  intersect 

3 3 


the  curved  lunar  surface  and  be  deflected  by  a second  scattering.  Further 
work  on  this  "bulk  shadowing"  effect  is  needed  to  bring  the  theory  of 
scattering  at  low  grazing  angles  into  better  agreement  with  the  experi- 
mental measurements  given  in  Figures  9 and  10.  It  is  worthwhile  to  point 
out  that  the  radius  of  curvature  of  the  Earth  is  nearly  four  times  that 

of  the  Moon,  so  the  "bulk  shadowing"  effect  would  be  of  reduced  importance 
for  the  Earth's  surface. 

From  this  brief  comparison  of  theory  and  experiment  several  ex- 
clusions can  be  drawn.  First,  the  quasispecular  scattering  theory  without 
shadowing  is  quite  acceptable  for  y > 15°  „ However,  for  y < 15°  shadowing 
Plays  an  increasingly  important  role  and  should  be  included  in  the  scatter- 
ing model.  For  15°  > y > 10°  the  shadowing  theory  of  Bass  and  Fuchs 
[Eqs.  (6)  and  (11)]  is  adequate  io  explain  the  experimental  data,  but 
for  y < loc  the  experimental  value  of  p2  falls  more  rapidly  than  the 
shadowing  theory  predicts.  This  failure  of  the  scattering  model  is 
thought  to  be  due  to  "bulk  scattering"  by  the  curved  surface  of  the  Moon. 

I uither  work  to  include  bulk  scattering"  in  the  model  of  Eq.  (k»0)  is 
needed  to  bring  theory  and  experiment  into  better  agreement  for  y < io° . 

Work  on  the  inclusion  of  "bulk  shadowing"  is  now  underway. 

To  sum  up  the  shadowing  theory  of  Bass  and  Fuchs  makes  a very 
necessary  improvement  in  the  scattering  model  at  low  grazing  angles 

(Y  < 15  however,  it  is  thought  that  "bulk  shadowing"  must  also  be 
included  for  grazing  angles  below  about  10°. 

D.  Relevance of  Scattering  Theory  to  Low-Angle  Radar  Tracking 

There  are  roughly  three  routes  by  which  improvements  in  scattering 
theory  will  manifest  themselves  as  improvements  in  low-angle  radar  per- 
formance. First,  an  accurate  scattering  theory  will  allow  accurate 
estimates  of  low-angle  elevation  errors  for  existing  radars,  since 


multipath  reflection  effects  are  the  major  source  of  error  at  low  angles 
Next,  a good  scattering  theory  will  enable  us  to  evaluate  new  schemes  for 
reducing  errors— for  example,  by  using  state  space  modeling  to  optimize 
the  radar  system  as  a whole.  Indeed,  features  revealed  by  a good  scatter- 
ing theory  may  well  suggest  new  methods  of  reducing  errors.  Finally,  a 
sufficiently  good  scattering  theory  might  enable  one  to  correct  at  least 
partially  for  multipath  errors  through  sophisticated  data  processing. 

1 . Error  Estimation 

lley  and  Parsons  (1955)  and  Evans  (1966)  make  elevation,-error 
estimates  using  specular  reflection  models  and  compare  their  estimates 
with  experimental  error  measurements.  While  the  theory  does  suggest  the 
main  features  of  the  measurements,  it  is  clearly  not  accurate  enough  to 
make  corrections.  Barton  (1974)  does  an  analysis  of  low-angle  elevation 
errors  using  a more  sophisticated  model  that  includes  "diffuse"  as  well 
as  specular  scattering.  He  concludes  that  approaches  exist  that  will 
allow  0.1  antenna  beamwidth  rms  accuracy  in  elevation  angle  at  elevation 
angles  as  low  a ' 4 beamwidth.  The  "diffuse"  scattering  theory  used  by 
Barton  (1974)  is  a quas ispecular  theory  modified  by  a "roughness  factor." 

An  experiment  has  been  proposed  by  Armstrong  et  al.  (1974)  to  tost  Barton's 
theory.  Obviously,  the  more  accurate  the  scattering  theory  is,  the  more 
accurate  will  be  the  error  estimates  for  a given  radar  system. 

2 , Evaluation  and  Synthesis  of  New  Techniques 

The  specular  reflection  theory  used  by  Hoy  and  Parsons  (1955) 
and  by  Evans  (1966)  suggested  to  them  techniques  for  error  reduction  that 
they  then  tried  in  practice.  For  example,  Evans  tried  erecting  12-ft-by- 
12-ft  shielding  screens  at  the  specular  point  to  reduce  multipath  errors. 
The  screens  did,  in  fact,  reduce  the  elevation  errors  by  about  a factor 
of  two  for  an  elevation  angle  of  about  2° . 
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Evans  (1966)  notes  that  during  experiments  at  X-band  frequencies 
the  mowing  oi  1-lt-long  grass  in  the  scattering  region  increased  the  mea- 
sured reflection  coefficient  by  a factor  of  about  4.  Such  interesting 
facts  have  suggested  that  modification  of  the  scattering  terrain  in  front 
oi  a radar  might  be  worthwhile  in  terms  of  error  reduction — especially 
something  as  easy  as  not  mowing  grass.  This  topic  is  discussed  at  more 
length  in  Section  111.  In  any  case,  a critical  theoretical  analysis  of 
*adai  siting  and  terrain  modification  as  useful  error— reduct ion  techniques 
requires  an  accurate  scattering  theory. 

A good  scattering  theory  is  also  required  for  an  overall  study 
of  a given  radar  system — for  example,  by  the  Kalman  filtering  approach 
discussed  in  Section  111-3-2-i  below. 

3 . Error  Correction 

A straightforward  acceptance  of  the  simple  specular  reflection 
model  used  by  Hey  and  Parsons  (1955)  and  Evans  (1966)  leads  one  to  expect 
that  a correction  algorithm  could  be  developed  for  multipath  errors.  Un- 
fortunately, the  situation  is  sufficiently  complicated  that  Evans  concluded 
such  a scheme  to  be,  in  general,  impractical.  More  recently,  White  (1974) 
has  successfully  developed  a system  that  greatly  reduces  the  effects  of 
multipath  over  water.  His  system,  which  models  the  multipath  effect  by 
assuming  specular  reflection,  is  discussed  in  more  detail  in  Section  III- 
B-2 • However,  were  a more  accurate  scattering  theory  available  it  would 
presumably  be  possible  to  use  the  theory  to  make  better  corrections,  or 
simply  to  use  the  incoming  data  in  a new  and. different  way.  Present  and 
future  electronic  hardware  will  make  possible  swift,  economical,  and  re- 
liable computation.  So  even  a rather  complicated  scattering  theory  could 
be  exploited. 
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III  RADAR  SYSTEMS  FOR  LOW-ANGLE  TRACKING 


A.  Gc  ne  ra 1 

There  are  many  system  configurations  and  system  techniques  that  may 
be  applied  so  as  to  achieve  good  low-angle  tracking  by  radar  systems. 
Depending  on  the  particular  application  and  the  desired  performance  level, 
various  combinations  of  these  techniques  may  be  employed.  In  Section  III- 
B-l  below,  the  general  radar  design  parameters  arc  examined  relative  to 
tracking  low-altitude  targets.  Section  III-B-2  examines  various  system 

techniques  that  may  be  applied  to  improve  both  new  designs  and  existing 
installations. 

In  much  of  what  follows,  it  is  assumed  that  a sufficiently  good  model 
exists  for  the  scattering  of  radar  signals  from  the  surrounding  terrain, 
lor  the  present  we  will  use  the  quasispecular  model  with  shadowing  dis- 
cussed in  Section  II  above.  However,  as  indicated  in  Section  II,  further 
work  is  needed  in  developing  a low-angle  scattering  model.  Such  develop- 
ments may  be  needed  before  some  of  the  techniques  in  this  section  can  be 
applied. 

B"  Kaclar  System  Design  for  Low-Altitude  Tracking  ,■  ility 

There  are  two  basic  approaches  that  can  be  applied  to  achieve  low- 
altitude  tracking  ability.  The  first  is  basically  a "brute  force"  approach 
in  which  the  basic  radar  is  designed  with  sufficient  bandwidth,  nower,  and 
antenna  aperture  to  separately  resolve  the  low-altitude  target  from  any 
possible  image  created  by  the  scattering  of  the  return  signal  from  the 
surrounding  surface.  The  second  approach  is  to  accept  the  fact  that  multi- 
path  scattering  of  the  target  return  signal  is  going  to  occur.  A model  of 
this  scattering  process  is  then  developed  and  used  in  the  processing  of 


43 


the  returned  signals  to  adjust  the  radar,  estimate  the  target  parameters, 
and  estimate  the  model  parameters  so  that  the  multipath  effect  can  be 
taken  into  account.  This  second  approach  may  be  applied  after  the  first 
approach  has  been  employed , or  it  may  be  applied  in  order  to  update  some 
existing  installations.  In  a recent,  useful  review  paper,  Barton  (1974) 
discusses  most  of  the  possibilities  of  the  first  approach  but  does  not 
cover  many  of  the  possibilities  of  the  second  approach.  Thus  che  emphasis 
of  this  chapter  will  be  primarily  on  the  second  approach.  One  should  bear 
in  mind,  while  considering  these  signal-processing  schemes,  that  large- 
scale  integrated-circuit  (LSI)  technology  has  made  and  will  continue  to 
make  (for  at  least  the  next  5 to  10  years)  data  processing  equipment 
smaller,  lighter,  cheaper,  and  more  reliable.  Thus  rather  sophisticated 
data-processing  schemes  that  might  have  been  prohibitive  in  the  past  due 
to  size,  expense,  or  reliability  considerations  may  well  be  possible  now. 

1 . Brute-Force  Methods  ■ 
a . Antenna  Aperture 

Perhaps  the  most  obvious  way  to  separate  a target  signal 
and  its  images  caused  by  scattering  from  the  terrain  is  to  increase  the 
antenna  aperture  Until  sufficient  angular  resolution  is  obtained  to  track 
the  target  and  reject  signals  scattered  by  the  terrain.  The  expression 
for  the  half-power  beamwidth  of  an  antenna  array  of  aperture  (length)  l 
is  given  by  Jasik  (1961,  p.  2-24)  as 

_ AX 

HP  Si  (24) 

where  A is  the  wavelength,  and  A is  the  constant  that  varies  slightly 

depending  on  the*  desired  sidelobe  level.  For  our  purposes,  A can  be 

chosen  as  1 where  8 is  in  radians.  This  corresponds  to  -25  dB  sidelobe 

HP 

levels  for  a Dolph-Tchebysche f f array.  If  f is  the  radar  frequency,  then 
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X = c/f 


(25) 


nnd 


c/f  i 


(26) 


N°VV  &HP  1S  a^so  aPProxiinately  the  angular  distance  from  the  antenna  pattern 

maximum  to  the  first  null;  so  setting  0 in  Figure  11  equal  to  0 we  have 

HP 


as  the  minimum  antenna  aperture  necessary  to  suppress  the  reflected  signal 
while  the  target  is  in  the  center  of  the  main  lobe.  It  is  theoretically 
possible  to  obtain  a considerably  narrower  main  lobe;  however,  all  such 
super  gain  designs  have  failed  in  practice  for  what  are  now  well  under- 
stood reasons,  as  explained  by  Buck  and  Gustincic  (1967).  The  angle  0 is 
easily  found  from  the  flat  Earth  radar/target  geometry  of  Figure  11  as 
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(28) 


The  present  assumption  about  the  scattering  process  is 

that  the  reflected  angle  y is  equal  to  the  incident  angle  y . Thus  we 
_i  3 i 

may  find  tan  *(h  /X  ) as  follows: 
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or 


tan  [ (h  - h )/G] 
t r 


The  antenna  aperture  i necessary  to  reject  the  reflected  return  signal 
may  now  be  determined.  The  results  of  this  calculation  are  illustrated 
in  the  last  column  of  Tables  1(a)  through  1(f)  and  were  produced  with 
the  computer  program  given  in  the  appendix. 

The  difficulty  of  constructing  large  antennas  grows  roughly 
3 

as  JL  (since  antennas  are  three-dimensional  structures).  Thus,  large  an- 
tennas become  very  expensive  and  impractical  for  very  large  apertures, 
lor  example,  from  Table  1 note  that  for  a range  of  54  km,  an  8-m  antenna 
aperture  is  required  for  a 10-GHz  radar  system  elevated  10  m,  tracking  a 
target  at  an  elevation  of  100  m.  Other  configurations  of  target  and  an- 
tenna height  are  illustrated  and  may  be  used  to  determine  the  feasibility 
of  improving  low-altitude  tracking  by  increasing  the  antenna  aperture. 


FIGURE  11  RADAR  ANTENNA/TARGET  GEOMETRY 


Table  1 


LOW-ALTITUDE  RADAR  SCATTERING  PARAMETERS 


(a)  Radar  frequency  (f)  = 10  GHz. 
Unidirectional  rms  slope  of 

scattering  surface  (tan  (30)  = 0.0, 
Radar  antenna  height  (hr)  = 10  m. 
Target  height  (lu)  = 100  m. 


Target 
' % 
Range,  r 

(m) 

Di  f ierentia 1 
Time  Delay, 6 

(ns) 

Range  to 
Sea  ttering 
Point,  r^ 

(m) 

Angle,  6 
( rad) 

Amplitude 

Ratio 

(Scattered/ 
Direct ) 

Required 
Antenna 
Aperture , J l 

(m) 

2,000 

3.3313 

181 

0.099914 

0.300 

4 ,000 

1.6672 

363 

0.049989 

KRgfK 

0.  599 

6,000 

1.  1116 

545 

0.033330 

KB 

0.  899 

8,000 

0.8338 

727 

0.024999 

0.925 

1.199 

10,000 

0.6670 

900 

0.019999 

0.940 

1.499 

12,000 

0.5559 

1090 

0.016666 

0.949 

1.798 

14,000 

0.4764 

1272 

0.014285 

0.956 

2.098 

16,000 

0.4169 

1454 

0.012500 

0.962 

2.398 

18,000 

0.3706 

1636 

0.011111 

0.966 

2.698 

20,000 

0.3335 

1818 

0.010000 

0.969 

2.998 

22,000 

0.3032 

2000 

0,009091 

0.972 

3.297 

24,000 

0.2779 

2181 

0.008333 

0.974  ' 

3,597 

26,000 

0.2565 

2363 

0.007692 

0.976 

3.897 

28,000 

0.2382 

2545 

0.007143 

0.978 

4.197 

10,000 

0.2223 

2727 

0.006667 

0.979 

4.497 

32,000 

0.2084 

2909 

0.006250 

n 980 

4.796 

34,000 

0.  1962 

3090 

0.005882 

0,981 

5.096 

36,000 

0.  1853 

3272 

0.005556 

0.  9§2 

5.396 

38,000 

0.  1755 

3454 

0.005263 

d.0§3 

5.696 

40,000 

0. 1667 

3636 

0.005000 

0.984 

5.996 

42,000 

0. 1588 

3818 

0.004762 

0.985 

6.295 

44,000 

0. 1516 

4000 

0.004545 

0.985 

6.595 

46,000 

0. 1450 

4181 

0.004348 

0,  986 

6.895 

48,000 

0. 1389 

4363 

0.004167 

0.  987 

7.195 

50 , 000 

0.1334 

4545 

0.004000 

0.987 

7,495 

52 , 000 

0. 1282 

4727 

0.003846 

0.988 

7.794 

54,000* 

0.1235 

4909 

0.003704 

0.988 

. E. 094 

* 

The  table  has  been  terminated  when  the  target  would  be  below  the  horizon 
for  a curved  Earth  with  standard  atmospheric  refraction. 


Table  i (Continued) 


LOW- Al.T ITVDE  RADAR  SCATTERING  PARAMETER: 


0>)  Radar  I ri  quenc\  II)  - 10  GH?. 

In  id  i ice  t i ona  1 .-ms  slope  of 

scattering  surface  (tan  • ) = 0 
Radar  antenna  height  (h,,)  = 10  m. 
Target  he  ight  (h.  ) 200  !t'. 


Targe ( 

( tinge,  i 

Ui  i t o n r.  1 1 a 1 
T i mo  Do  1 a y , 

ftungo  to 
.^ca  t tor  i 14* 
Point , r j 

Angle,  - 

\mpl i tude 
Ha  t i o 

Required 

Antenna 

Aperture 

(m) 

(ns) 

{ m) 

(rad) 

(Sea , to  red  / 
l>  recti 

(m) 

1 

1 , 

6 , 

8, 

1 0 , 000 
12,000 

1 1 . 000 
1(5,000 

1 8 . 000 
20 , 000 
22,000 
2 A ,000 
26,000 
28 , 000 
.10,000 

32.000 
34 , 000 

36.000 

38.000 
-10,000 
-12,000 
-hi  ,000 
-16,000 
-18 , 000 
DO , 000 

52 . 000 

5- 1 , 000 

56.000 

58.000 
(50,000 

62.000 

6- 1,000 
66 , 000 
68,000 
70.000" 


6,6379 
3.33  5 
2 , 2 -2-1 
1 . 6672 
9 

1 , 1 1 1 (5 
0,  9529 
0. 8338 
0.7-111 
0. 6670 
0.6061 
0. 5559 
0. 513) 
0.4761 
0,-1117 
0.  -1169 
0.392-1 
0. 3706 
0,351  1 
0.3335 
0, 3176 
0,3032 
0.2900 
0, 2 779 
0,2668 
0,2565 
0,2470 
0,  2,382 
0,2300 
0,2223 
0,2151 
0.2084 
0.2021 
0 , I 962 
0,  1906 


0.199332 
0.099916 
0.066642 
0. 049990 
0.039995 
0.033330 
0. 028569 
0.  02 1990 
0,022221 
0.019999 
0. 018181 
0, 010666 
0. 015384 
0,014285 
0.01 333.3 
0,012500 
0.01 1765 
0,01 1111 
0.010526 
0.010000 
0,009524 
0,009091 
0.008696 
0.008333 
0. 008000 
0,007692 
0,007107 
0,007143 
0.006897 
0. 006667 
0. 006452 
0.006250 
0 , 00606  I 
0.  005882 
0.005714 


O.  575 
0.  751 
0,824 
0.864 
0,889 
0 , 906 
0 9)9 
0,929 
0 , 936 
0.9-12 
0 94  7 
0.951 
0 , 955 
0.958 
0 . 96 1 
0.  963 
( .965 
0.967 
(, . 9 


0.  150 
0.  300 
0.-149 
0.  599 
0.749 

0.  899 
,019 

1.199 

1 .  : 9 
1,499 
1.648 
1. 798 
1 .94  8 
2,098 
2,248 
2,398 
2.548 
2.698 
2. 848 
2.998 
3.1-17 
3.297 
3.447 
3 . 597 
3.74  7 
3.897 
4.04  7 
4.197 
4.34  7 
4.  197 
4.6  16 
4.796 
4.946 
5 . 096 
5.246 


111'.-  table  has  been  term! nu  led  Alien  the  target  would  be  below-  the  lioil/.on 
lor  a curved  Earth  with  standard  atmospheric  re  1 rue  t i on . 
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T.'i  h le  1 (Continued) 


1.0VV-A LT I Tl  UK  RADAR  SCATTKH1NG  PARAMETERS 


(c)  Radar  I requoncy  (I)  10  Ollz. 

Rnldi  rod  lonn  1 rms  slope  ol 

sen  t to  ling;  surface  (tan  f!r))  - 0.0. 
Radar  antenna  height  (h  ) = 10  m. 
Target  height  (h()  300  m 


T n rge  t 
Ra  nge , r 

(ra) 

LM  f fe rent  i a i 
T i me  Delay. 

(ns) 

Hanin>  tn 
Sent  tel  l llg 
.’iiint  , r ! 

(in) 

Anglo,  9 
(rad) 

Ampl 1 1 udo 
Ra  1 i o 

(Seal  to  rod/ 
Direct ) 

Retfu  1 rod 
Antenna 
Aperture,  l 

(m) 

2.000 

9.8958 

64 

0.297773 

0. 455 

0 . 1 00 

1 , 000 

1.9893 

129 

0. 119719 

0.660 

0.200 

6,000 

3,3310 

193 

0.099917 

0.751 

0.300 

H , 000 

2,4999 

258 

0.07  1965 

0.  807 

0.399 

10,000 

2 . 0004 

322 

0.059982 

0.842 

0.499 

12,000 

1 ,6672 

387 

0.049990 

0.866 

0.599 

14 ,000 

1 .1291 

451 

0.042851 

0.  88.8 

0.699 

16,000 

1.2506 

516 

0.037496 

0.  897 

0.799 

18,000 

1,1116 

580 

0.033330 

0.  908 

0.899 

20,000 

1 . 0005 

645 

0.029998 

0.916 

0 . 999 

22 , 000 

0.9096 

709 

0.027271 

0.  924 

1 . 099 

2-1,000 

0.8338 

771 

0.024999 

0.930 

1 . 199 

26 . 000 

0,7696 

838 

0.023076 

0.935 

1.299 

2 8 , 000 

0.7147 

903 

0.021428 

0.939 

1 . 399 

30,000 

0.6670 

967 

0.019999 

0.943 

1 . 499 

32 , 000 

0.6253 

1032 

0.018749 

0.947 

1 . 593 

31,000 

0,5886 

1096 

0.0176-17 

0.950 

1.698 

36 , 000 

0.5559 

1 161 

0.016666 

0.952 

1,798 

38 , 000 

0,5266 

122  5 

0.015789 

0.955 

1 .898 

40,000 

0 . 5003 

129 

0.015000 

0.957 

1.998 

42  000 

0.4764 

i:  si 

0.014285 

0.959 

2.098  ’ 

44,000 

0.4548 

1 1 19 

0.013636 

0.961 

2.198 

16,000 

0.4350 

1 483 

0.013043 

0.962 

2.298 

-18,000 

0.4169 

1518 

0. 0i2500 

0.964 

2.398 

50 , 000 

0.4002 

1612 

0.012000 

0.965 

2,498 

52 , 000 

0.3848 

1677 

0.011538 

0.  966 

2 . 598 

54 , 000 

0.3706 

1741 

0,011111 

0.968 

2.693 

56 , 000 

0.3573 

1806 

0.010714 

0.969 

'.798 

58 , 000 

0.3450 

1870 

0.010345 

0.970 

2.898 

60,000 

0.3335 

1935 

0.010000 

0.97  1 

2.998 

62,000 

0. 3227 

2000 

0.009677 

0.  972 

3.098 

61,000 

0. 3127 

20C4 

0.009375 

0.973 

3,197 

66,000 

0.  3032 

2129 

0.009091 

0.973 

3.297 

68,000 

0.2943 

2193 

0.008823 

0.974 

3.897 

70,000 

0.2859 

2258 

0.008571 

0.9,5 

8.497 

72,000 

0.2779 

2322 

0.008333 

0.976 

3.597 

74,000 

0.2704 

2 387 

0.008108 

0.976 

3.697 

76 , 000 

0,2633 

2 151 

C .007895 

0.977 

3.797 

78,000 

0,2565 

2516 

0.007692 

0.977 

3,897 

80,000 

0.2501 

2580 

0.007500 

0.  978 

3.997 

82,000 

012440 

2645 

0.007317 

0.978 

4.097 

84,000* 

0.2382 

2709 

0.007143 

0.979 

4.197 

Tliu  table  has  been  tc‘i*minated  when  the  target  would  be  below  the  horizon 


lor  a curved  Earth  with  standard  atmospheric  refraction. 


Table  1 (Continued) 


LOW-ALTITUDE  RADAR  SCATTERING  PARAMETERS 

(d)  Radar  freqi  (f)  - 10  GHz. 

Unidirectio  ;cl  .ms  slope  of 

scattering  surface  (tan  0Q)  = 0.05. 
Radar  antenna  height  (hr)  = 20  m. 
Target  height  (ht)  = 100  m. 


Ta rget 
Range , r 

(m) 

Dif  ferent  ia  1 
Time  Delay,  5 

(ns) 

Range  to 
Scattering 
Point,  l-j 

(m) 

Angle,  9 
(rad) 

Amplitude 

Ratio 

( Scattered/ 
Direct ) 

Required 
Antenna 
Aperture,  i 

(m) 

20,000 

0.6671 

3,333 

0.010000 



0.102 



2 . 998 

40,000 

* 

0.3335 

6 , 666 

0.005000 

0.052 

5.996 

o 

o 

o 

o 

3 

0.2223 

10,000 

0.003333 

0.035 

8.994 

(e)  Radar  frequency  (f)  = io  GHz. 
Unidirectional  rms  slope  of 

scattering  surface  (tan  pQ)  = 0.05. 
Radar  antenna  height  (hr)  = 60  m. 
Target  height  (ht)  = 300  m. 


T a rge  t 
Range,  r 

(m) 


20,000 

40.000 

60.000 
80,000* 


Differential 
Time  Delay,  6 

(ns) 


2.0011 
1.0006 
0.6671 
0. 5003 


Range  to 
Scattering 
Point,  r^ 

(m) 


1,250 

2,500 

3,750 

5,000 


Angle , 8 


(rad) 


0.029998 
0.015000 
0.01 0000 
0.007500 


Amplitude 

Ratio 

(Scattered/ 

Direct) 


0.252 

0.134 

0.091 

0.069 


Required 
Antenna 
Aperture,  l 

(m) 


0.999 

1.998 

2.998 
3.997 


The  table  has  been  terminated  when  the  target  would  be  below  the  horizon 
ior  a curved  Earth  with  standard  atmospheric  refraction. 
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Table  1 (Concluded) 


LOW-ALTITUDE  RADAR  SCATTERING  PARAMETERS 


(f  ) Radar  frequency  (f)  = 10  GHz. 
Unidirectional  rms  slope  of 

scattering  surface  (tan  0O)  = 0.05. 
Radar  antenna  height  (hr)  = 60  m. 
Target  height  (ht)  = 300  m. 


Target 
Range , r 

(m) 

Differential 
tme  Delay,  6 

(ns) 

Range  to 
Scattering 
Point,  r^ 
(m) 

Angle,  0 
( rad) 

Amplitude 

Ratio 

(Scattered/ 

Direct) 

Required 
Antenna 
Aperture,  l 
(m) 

20,000 

6.0033 

3,333 

0.029997 

0.279 

0.999 

40,000 

3.0019 

6 , 666 

0.015000 

0. 150 

1.998 

60,000 

2.0013 

10,000 

0.010000 

0. 102 

2.998 

80,000 

1.5009 

13,333 

0.007500 

0.077 

3.997 

100,000 

1.2007 

16,666 

0.006000 

0.  062 

4.996 

The  table  has  been  terminated  when  the  target  would  be  below  the  horizon 
for  a curved  Earth  with  standard  atmospheric  refraction. 


One  way  of  increasing  the  effective  antenna  aperture  with- 
out actually  constructing  a larg^  dish  antenna  is  to  mount  a series  of 
small  antennas  on  a receiving  tower.  This  method  is  similar  to  the  aper- 
ture-synthesis technique  familial'  in  radio  astronomy — cf.  Christiansen 
and  Hogbom  (1969),  or  Kraus  (1969),  Such  a receiving  tower  would  be  ef- 
fective only  if  the  echo  signal  were  sufficiently  strong  and  improved 
resolution  were  required  only  in  elevation  (to  discriminate  against  the 
image  echo  signal  of  Figure  1),  In  its  simplest  realization,  suc'i  a tower 
would  be  about  100  m high,  and  would  have  several  (one  to  ten)  receivers 
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r 


\ 


«"»*  vertically  along  the  tower.  The  tower 
""iC"  ”y  " "ay  “0t  b«  ‘“-‘“d  near  the  transmitter,  would  not  do  dim-  ' 
cult  to  construct  and  put  in  place,  the  main  engineering  requirement  would 
position  of  each  of  the  receivers  would  have  to  be  known  to  a 
fraction  oi  the  wavelength.  Of  the  tower  swayed,  a correction  could  he 
»nde  electronically  11  the  .agniu.de  oi  the  swaying  were  .ensured  in  real 
«i»c.)  .Such  a receiver  could  give  a resolution  in  the  vertical  direction 
of  about  X/d  = 10  rad  for  a 100-.  lower  at  !0  The re  wou 1 d , of 

course,  he  no  useful  resolution  in  the  horizontal  direction.  For  i.provod 

th°  P"1CO  °f  * increase  in  co.pJoxity,  the  antennas 

could  have  narrower  boa.s  and  ho  .ado  to  sweep  in  the  horizontal  direction 

in  synchronization  with  the  transmitter  nr.  + 

transmitter  Or,  the  tower  could  be  made 

higher  with  more  receiving  antennas. 

It  is  also  possible  to  ny  tho  supplomentary  recelver  from 
a balloon,  rocket,  airplane,  or  other  airborne  support,  ngain  only  if  its 
position  is  constantly  monitored.  Such  monitoring  could  probably  be 
achieved  with  three  narrowband  microwave  transponders,  chosen  so  as  not 
to  interfere  with  tho  radar  return,  and  a simple  computer.  If  , rockot 
wore  used,  It  would  he  launched  only  upon  receipt  o,  an  in, cresting  signal 
With  a receiver  that  would  slowly  drift  down  on  a parachute  or  balloon 

dat"  C0“Id  b°  ,0l0m"0rM'  - ™ . »‘w.  Tho  main  disadvantage 

Of  a si„g,o  receiver  (or  of  a small  number  of  receivers)  is  that  the  high 

elution  in  the  vertical  direction  would  suffer  from  high  sidelobes, 
which  might  require  sophisticated  processing.  „ ,«r,  mapy 

co. vers,  these  sidelobes  could  be  made  negligibly  small  through  appropriate 
apodiza  tion. 
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b. 


Radar  Operation  at  mm  Wavelengths 


In  Eq.  (24)  we  note  that  the  half-power  beamwidth  (9)  of 
an  antenna  is  given  by  0 % \/i.  The  previous  section  noted  the  improvements 
in  low-angle  radar  performance  that  could  be  made  by  reducing  0 through  an 
increase  in  the  antenna  aperture  (length)  JL.  Similar  improvements  could 
also  be  made  by  reducing  the  radar  operating  wavelength  X,  thus  avoiding 
the  difficulties  of  large  antenna  structures.  However,  at  frequorcies 
above  40  GHz  atmospheric  attenuation  begins  to  play  an  important  role, 
and  transmitting  and  receiving  systems  become  less  efficient.  This  and 
other  applications  of  mm  wavelength  devices  have  been  investigated  by 
R.  0.  Hundley  (1975). 

Evans  and  Hagfors  (1968,  p.  391)  give  a convenient  expres- 

sion  for  the  expected  signal-to-noise  ratio  of  a radar  system.  We  can 

easily  modify  this  expression  to  suit  our  purposes  by  introducing  the 

factor  exp(-2o<r)  to  account  for  atmospheric  absorption.  Letting  the 

2 

antenna  gaii  , G = 4rcA  A , we  have 

e 


where 


2 -2orr 
P P A ct  e 
r _ t e 

p 4 2 

N 4rcr  X kT  B 
s 


~ received  echo  power 


P = receiver  system  noise  oower 
N 

P = peak  transmitted  power 


A^  = effective  antenna  area  = n(Z/2)‘ 


° = target  cross  section 


o = atmospheric  attenuation  coefficient 
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r = range  to  target 


A.  = radar  system  wavelength 

= receiver  system  noise  temperature 
B = receiver  system  bandwidth 

and  all  are  in  SI  (MKS  units).  Suppose  we  now  require  that  our  system 

have  a 10  dB  signal-to-noise  ratio  for  low  angle  tracking,  i.e. , (P  /P  ) = 10 

r N 

Knowing  the  other  parameters  we  can  solve  for  the  required  transmitter  peak 
power,  P^_, 

Taking  nominal  system  parameters  as  follows: 

2 

Ae  = = JV/4,  i.e.,  antenna  diameter,  i = 1 m. 

2 

0-  = 0,  1 m 

B =0.1  GHz 

and  using  the  above  equation  we  can  construct  Table  2 below.  The  factors 
governing  the  parameter  choices  are  discussed  in  due  course. 


Table  2 

LOW-ANGLE  RADAR  TRACKING  AT  mm  WAVELENGTHS 


Ta  rge  t 
Ra  nge , r 

Radar 

Wavelength,  \ 

Atmospheric 
At  tenuation 
Coefficient,  a 

System 

Noise 

Temperature,  T 

s 

Antenna 
Beamwidth,  9 

Required 
Transmitter 
Power,  P 

t 

(m) 

(mm) 

(m  '*’) 

(K) 

(rad) 

(W) 

20,000 

10mm  (=30GHz) 

-5 

1.6  X 10 

650 

0.01 

4 

4.4  X 10 

50,000 

10mm  (=30GHz) 

1.6  X io  5 

650 

0.01 

6 

4.4  X 10 

20,000 

1.35mm  (=220GHz) 

-4 

2.1  X 10 

2,700 

0.00135 

6 

7.7  X 10 

50,000 

1.35mm  (=220GHz) 

2.1  X 10~4 

2,700 

0.00135 

9.0  X 1013 
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Using  Table  1 we  can  compare  the  low-angle  tracking  per- 
formance of  the  K band  (30  GHz)  and  mm  band  (220  GHz)  radars  considered 
a 

in  Table  2.  Looking  in  Table  1(a),  for  example,  we  find  that  for  a radar 

situated  10  m above  a flat  Earth  and  tracking  a target  flying  at  100  m 

altitude  the  K band  radar  of  Table  2 will  track  out  to  a range  of  about 
a 

20  km.  The  mm  wave  radar  by  contrast  will  track  virtually  out  to  the 
horizon  some  50  km  away.  This  difference  in  performance  is,  of  course, 
due  to  the  extremely  narrow  beam  pattern  of  our  1-m  diameter  dish  at  mm 
wavelengths.  The  superior  performance  of  the  mm  wave  radar  is  not 
achieved  without  effort  however.  The  peak  power  required  (P  ) is  ~ 200 
times  greater  for  the  mm  wave  radar  at  r = 20  km  and  increases  dramati- 
cally as  r increases.  The  larger  values  of  P^  required  at  the  mm  wave- 
length are  due  to  increased  receiver  noise  (higher  T ) and  more  importantly 
to  the  much  greater  amounts  of  atmospheric  absorption  (higher  a).  Atmos- 
pheric absorption  varies  strongly  with  frequency  and  the  1.35-mm  wavelength 
is  near  a relative  minimum.  The  r = 20  km  values  of  P^  required  in  Table  2 
are  commensurate  with  current  K band  technology;  but  peak  powers  reliably 
available  at  around  1-mm  wavelength  fall  in  the  1 to  10  watt  region,  far 
below  the  requirements.  We  can  therefore  conclude  that  a significant 
improvement  in  low-angle  radar  performance  can  be  obtained  by  going  to 
a 10-mm  wavelength  (K  band)  within  current  technological  limits.  However, 
further  significant  improvements  obtainable  by  going  to  ~ 1-mm  wavelengths 
will  require  a very  large  improvement  in  transmitter  peak  power.'  One  can 

use  pulse  compression  techniques  to  reduce  the  required  peak  power  by  a 
4 

factor  of  up  to  10  . Even  using  pulse  compression,  substantial  improve- 
ments in  transmitter  peak  power  must  be  realized  to  exploit  the  superior 
low-angle  performance  of  a ~ 1-mm  wavelength  radar. 


The  numbers  ciuoted  in  Table  ? and  the  conclusions  above 
are  dependent  on  the  parameter  choices  made,  and,  in  particular,  the 
values  of  a and  . Atmospheric  absorption  is  due  mainly  to  water  vapor 
and  oxygen  molecules.  The  values  of  O'  used  are  for  1%  water  vapor  and 
an  atmospheric  pressure  of  760  mm  of  Hg.  Absorption  peaks  due  to  molecu- 
lar resonance  bands  make  a a strong  function  of  frequency.  The  particular 
operating  wavelengths  in  Table  2 correspond  to  relative  minima  in  atmos- 
pheric absorption.  Skolnik  (1962,  p.  516)  gives  the  plot  of  a versus 
wavelength  from  which  the  values  of  Table  2 were  taken. 

The  system  noise  temperatures  used  in  Table  2 are  based  on 

the  use  of  microwave  integrated  circuits  (MICs).  While  current  MICs  at 

30  GHz  (10  mm)  can  achieve  a 5-dB  noise  figure  (F) , no  MICs  at  220  GHz 

(1.35  mm)  have  yet  been  constructed  to  the  authors'  knowledge.  At  the 

higher  frequency,  workers  in  the  field  think  a 10-dB  noise  figure  is 

achievable  and  lower  values  could  probably  be  achieved  with  sufficient 

effort.  The  system  noise  temperature  T was  calculated  using  T = T (F-l) 
o s so 

where  T°  = 300  K.  As  can  be  seen  from  the  equation  above,  the  peak  power 

required  (P^  is  directly  proportional  to  Tg;  so  the  realization  of  a 

x . 35-mm  radar  could  well  be  aided  by  lowering  the  system  noise  temperature 

through  low-noise  receiver  development. 

Another  factor  that  will  significantly  influence  mm  wave 
radar  performance  is  the  accuracy  to  which  the  radar  antenna  can  be 
constructed.  At  short  wavelengths  antenna  surface  tolerances  become 
very  close.  Ruze  (1952)  considered  the  effects  of  parabolic  dish  dis- 
tortion on  antenna  gain.  For  a worst  case,  the  gain  reduction  is  given 
by 
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who  re 


a = ±2* 

\ 

and  d = rms  surface  deviation  from  true  parabolic  shape.  Since  the 
antenna  gain  (G)  is  related  to  the  antenna  beamwidth  (0)  by  G = 4tc02, 
any  significant  reduction  in  antenna  gain  will  significantly  reduce 
low-angle  tracking  performance.  For  example,  suppose  we  can  construct 
an  antenna  surface  to  a tolerance  (d)  of  0. 1 mm  (~  4 X 10~3  in.).  At 
X = 1°  mm  (30  GHz>.  thc  Sain  reduction  is  about  0.1  dB  and  quite  ac- 
ceptable. If  we  try  to  use  the  same  antenna  at  \ = 1.35  mm  (220  GHz), 
the  gain  reduction  is  about  9 dB  and  8 increases  by  about  3 times, 
which  could  well  be  important. 

In  summary,  mm  wavelength  radars  can  achieve  good  low- 
angle  tracking  performance  by  allowing  one  to  obtain  very  narrow  antenna 
beamwidths  with  relatively  small  antennas,  such  as  the  1-m  diameter  dish 
considered  in  Table  2.  However,  as  one  operates  at  progressively  higher 
frequencies,  two  factors  require  that  the  peak  transmitter  power  P be 
progressively  larger.  These  factors  are  the  higher  atmospheric  attenua- 
tion and  higher  receiver  noise  temperature  at  higher  frequencies.  While 
the  required  effective  peak  power  at  10  mm  is  within  the  current  tech- 
nology, Pt  rises  well  above  a megawatt  at  1.35  mm  for  a target  at  20-km 
range.  Such  a power  level  requires  substantial  increases  in  transmitter 
peak  power  capability  at  wavelengths  around  1 mm  in  order  to  exploit  the 
very  much  improved  low-angle  tracking  performance  achievable  at  these 
wavelengths.  Improvements  in  receiver  system  noise  temperature  over  the 
estimates  used  in  Table  2 would  also  be  helpful  in  reducing  the  required 
value  of  Pt,  perhaps  by  as  much  as  a factor  of  4.  It  should  also  be 
noted  that  at  short  wavelengths— e.g.  , 1.35  mm— antenna  surface  tolerances 
become  very  tight  indeed  (~  0.05  mm)  if  satisfactory  performance  is  to 
be  achieved. 
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Range  Resolution 


c . 

If  the  target  must  be  tracked  while  the  ; ntenna  pattern 
includes  both  the  target  and  the  horizon  (due  to  insufficient  antenna 
aperture),  then  the  ability  to  determine  the  target  altitude  may  be  a 
function  of  the  range  resolution  of  the  system.  This  occurs  sihep  range 
gates  may  be  used  in  the  tracking  circuits  to  exclude  the  reflected  signal 
if  it  can  be  sepa~„ted  in  time  from  the  direct  target  signal. 

Thus,  sufficiently  short  pulses  will  produce  separate  direct 
and  scattex'ed  returns.  Short  pulses  are  associated  with  large  bandwidths 
and  it  is  easily  shown  [by  Barton  (1974),  for  example]  that  low-altitude 
tracking  ability  is  proportional  to  the  effective  radar  system  bandwidth, 
whether  it  is  implemented  directly  as  suggested  above  or  by  pulse  compres- 
sion, frequency  agility,  FM  sweeps,  etc.  Large  system  bandwidths  are 
usually  expensive,  and  for  many  systems  of  interest  it  is  often  insuffi- 
cient to  provide  the  desired  low-altitude  tracking  capability. 

Tiie  difference  in  arrival  times  of  the  direct  and  the  re- 
flected return  is  determined  by  the  antenna,  ground,  and  target  geometry. 
The  difference  in  path  length  for  a flat  Earth  can  be  determined  from 
Figure  11,  as  shown  in  Eqs.  (33)  and  (34). 


The  direct  path  length  is 


r 


V°2  + 


(33) 


while  the  reflected  path  length  is 


(34) 
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Thus,  the  path  difference,  6,  is 


However,  G 
given  by 


X + T and  h /X  = h /T 
£ r s t 
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h /(h  + h ) 

r t r 
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or 
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(h  + h )2 
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(35) 


(36) 


and  the  time  difference  is  6/c  where  c is  the  velocity  of  light.  The 
above  function  has  been  calculated  and  is  given  in  the  second  column  of 
Tables  1(a)  through  1(f)  for  various  parameters  of  the  problem.  The  most 
a pparent  feature  of  the  delay  time  is  that  for  any  target  at  long  range, 
the  difference  in  arrival  times  is  very  small,  on  the  order  of  0.1  ns. 
Ihus,  oven  extremely-widc-bandwidth  radars  could  not  directly  separate 
the  two  signals  in  time. 


d . Siting 


The  last  brute-force  method  concerns  the  selection  of  a 
proper  site  for  the  radar.  From  Tables  1(a)  through  1(F)  it  is  apparent 
that  a sufficiently  large  antenna  height  will  provide  a->  desired  low-angle 
tracking  capability.  Very  tall  antennas  are  not  always  practical,  so 
this  technique  is  useful,  but  limited.  The  second  site  consideration  is 
the  surrounding  terrain.  Clearly,  a perfectly  smooth  terrain  would  be  a 
very  good  reflector  of  radar  signals,  while  suitably  rough  surfaces  could 
virtually  eliminate  the  specularly  scattered  signal.  Tables  1(d),  (o), 
and  (f)  also  illustrate  the  effect  of  surface  roughness  on  the  reflected 
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signal  amplitude  in  the  fifth  column.  The  scattering  model  used  is  that 
of  Eq.  (20)  with  C = D = 1 and  e =3  in  Eqs.  (22)  and  (23).  Note  that  a 
unidirectional  rms  surface  slope  tan  of  only  0.05  greatly  attenuates 
the  reflected  return  signal  [see  column  5 of  Table  1(a)  through  1(f)]. 

This  calculation  uses  a very  simple  model,  however,  and  as  extensively 
discussed  in  Section  II,  all  known  low-angle  scattering  models  are  suspect. 

2-  Compensation  Techniques  for  Low-Altitude  Tracking 

Once  the  basic  radar  parameters  of  frequency,  bandwidth,  antenna, 
aperture,  and  site  have  been  fixed,  then  additional  techniques  may  be  ap- 
plied to  compensate  for  the  reflected  return  signal. 

a • Signal  Processing  and  Measurements 

The  purpose  of  this  subsection  is  to  demonstrate  that  the 
difficulties  of  tracking  a low-altitude  target  are  mostly  an  artifact  of 
the  paiticular  data-processing  algorithm  a monopulse  radar  uses,  and  that, 
with  an  appropriate  but  modest  change  in  the  monopulse  radar  configuration 
and  processing  algorithm,  a radar  can  successfully  track  targets  at  low 
elevation  angles. 

Any  processing  algorithm  presumes  a definite  model  of  the 
reflection  coefficient  from  the  ground  (e.g.,  specular  reflection  and 
geometrical  optics),  and  it  is  important  to  assess  how  sensitive  such  an 
algorithm  is  to  deviations  from  the  assumed  reflection  model.  Such  devia- 
10ns  will  surely  occur  in  real  systems.  For  example,  even  smooth  mirror 
surfaces  (the  ideal  specular  reflection)  have  corrections  due  to  thysical 
optics — i.e.,  Fresnel  zones,  when  the  radar  is  situated  within  a rew  an- 
tenna diameters  of  the  surface. 
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The  key  question  in  the  discussion  to  follow  is,  how  many 
unknowns  there  are  in  the  reflected  radar  signal,  and  how  many  measure- 
ments the  radar  must  make  to  determine  the  unknowns.  Clearly,  in  a 
multipath  low-elevation-angle  situation,  the  number  of  unknowns  increases 
due  to  the  reflected  wave.  The  radar,  will  then  have  to  make  more  measure- 
ments to  obtain  information  about  the  reflecting  surface  to  identify  the 
true  elevation  of  the  target. 

In  general,  the  idea  is  to  express  the  properties  of  the 
reflecting  surface  according  to  a model  that  is  specified  by  a small  numbei 
of  adjustable  parameters  (usually  two  or  three).  The  numbe"  of  receive 
beams  on  the  radar  is  increased  to  permit  two  or  three  addit  onal  measure- 
ments, v.'hich  then  allows  one  to  determine  the  adjustable  perameters  in  the 

reflection  coefficient  ant1  hence  determine  the  bearing  of  the  target. 

First,  let  us  see  how  these  considerations  work  in  a stan- 
dard monopulse  tracking  scheme  (Figure  12).  Two  signals,  and  , are 
received  whose  amplitude  and  phase  are  given  by 

9 j 1/  p 

51  = A C(0T  " 9I1  ~ A9)  ° (37) 

52  = A G(6t  - 0B  + A0)  e2ik1’  (38) 

whore  r denotes  the  range;  A is  a complex  amplitude  determined  by  the  range 
the  phase  shift  on  reflection,  the  radar  cross  section  of  the  (argot,  and 
other  factors;  and  G is  the  gain  pattern  of  the  radar  antenna.  Thus,  there 
are  four  unknowns,  r,  Rc(A),  In.’(A),  and  6 , and  four  measurements,  the  real 
and  imaginary  parts  of  and  , which  allow  us  to  determine  the  four  un- 
knowns. In  particular,  one  employs  the  algorithm 
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to  find  the  error  signal  that  locates  the  target. 


FIGURE  12  RECEIVE  BEAMS  OF  A MONOPULSE  RADAR.  The  upper  beam  is  denoted  by  1, 
while  2 denotes  the  lower  beam.  0R  gives  the  boresight  elevation  angle  of  the  radar, 
and  the  upper  and  lower  beams  are  separated  from  the  boresight  direction  by  an 
angle  AO. 


Now  consider  a three-beam  system  and  reflection  from  the 

th 

ground  as  shown  in  Figure  13,  The  i signal  will  then  be 

S.  = A Gv0  - 9 ) e2lKr  + 

l 1 T i 

where  represents  the  radiation  from  the  target  directly  into  the  antenna 
and  A (0 ) the  radiation  that  is  received  after  ground  reflection.  The  re- 
flection coefficient  is  p ( 0 ) ei'lf(9)> 
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FIGURE  13  A THREE-BEAM  RADAR.  St  ch  a system  can  distinguish  between  a target  and  its 
reflected  image.  0B  is  the  br resight  angle,  and  0.  (i  = 1,  2,  3)  denote  the  pointing 
directions  of  the  radar  beams. 


It  is  evident  that  if  nothing  is  assumed  to  be  known  about 

the  reflection  coefficient,  then  each  new  measurement  S (really  two  mea- 

i 

surements  amplitude  and  phase),  corresponding  to  a new  value  of  8 , will 
only  give  two  new  pieces  of  information  about  the  function  A(8)p(8)  e1^0), 
which,  in  prim  pie,  has  an  infinite  number  of  degrees  of  freedom.  In  such 
circumstances,  we  can  never  measure  enough  about  the  reflection  coefficient 
to  determine  the  target  bearing. 


Consequently,  a reflection  coefficient  model  must  be  adopted. 
The  one  mainly  used  so  far  is  the  specular-ref lection/geometrical-optit 
approach  whereby  _ 


ICS 


. itii  (0 ) i>) 

p(8)  e.-  } = P o 6(6  + 9 - 20  ) 
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(41) 


where  p is  the  reflection  coefficient,  * is  the  phase  shift  of  the  reflected 
* 

ray  relative  t?  the  direct  ray,  6(  ) is  the  Dirac  delta  function,  and 
* 

9p  ls  thc  inclination  of  the  reflecting  plane  to  the  horizontal.  This 
equation  represents  the  reflection  coefficient  by  a three-parameter  model. 

In  over-water:  applications,  the  parameters  may  be  reduced  to  two  by  the 
assumption  8;  = 0. 
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Each  signal  S is  then  given  by 
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We  can  suppose  A^  and  A^  to  be  real  by  absorbing  any  phase  shifts  into  r 

and  t.  In  this  case,  we  have  six  unknowns,  r,  A , A p,  i|f,  6.  and  6 , and 

12  T p 

six  measurements  of  amplitude  and  phas  of  S. , i = 1,  2,  3.  In  principle, 
9 can  be  determined 

White  (1974;  see  also  Section  IV)  has  used  a restricted 
version  of  this  scheme  over  water  <9  = 0 is  assumed)  wherein  only  the 

real  parts  of  the  signal  are  used. 

The  three  unknowns  are  then  A^  cos  kr,  A p and  (kr  + \|r), 
and  9^,  and  the  three  measurements  are  Re(S.)  = 1,  2,  3.  Again,  9 can 
be  found.  White’s  scheme  shows  that 


Error  signal  a 


(43) 


when  9t  is  much  less  than  a beamwidth,  so  that  the  signal-to-noise  ratio 
is  low  when  the  elevation  angle  is  low.  But  it  is  very  important  to  point 
out  that,  this  processing  algorithm  keeps  the  errors  bounded  and  the  beam 
on  the  target. 

It  should  be  quite  straightforward  to  modify  existing 
monopulse  radars  to  have  low-altitude  tracking  capability.  Figure  14 
shows  how  feeds  can  be  rearranged  to  achieve  three  vertical  beams. 

We  have  always  supposed  the  reflection  coefficient  to  be 
large  (P  1 ) . In  a particular  air  Jefense  situation,  the  radars  may  be 
located  in  forest  or  other  areas  where  the  scattering  is  weak  and  diffuse 
leading  to  very  little  multipath.  Our  proposed  schemes  will  work  regard- 
less of  the  value  of  p, 
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STANDARD  LOW-ELEVATION 

MONOPULSE  TRACKER 


1 ► HORIZONTAL 

FIGURE  14  HORN  ARRANGEMENT  ON  STANDARD  MONOPULSE,  AND  LOW-ELEVATION 
TRACKER.  The  numbers  refer  to  the  antenna  beams  for  the  standard  monopulse 
(Fig.  12)  and  three-beam  (Fig.  13)  radars. 

Tn  sunmary,  modest  modifications  of  existing  radars  will 
give  them  low-altitude  tracking  capability.  But  more  work  is  required 
to  find  the  optimum  algorithms  and  to  determine  their  sensitivity  to 
various  reflection-coefficient  models.  The  further  additional  advantages 
of  radar  netting  are  . .ined  below. 

b . Radar  Nets  and  Bistatic  Operation 

The  netting  together  of  several  radar  installations  is  a 
useful  technique  for  improving  overall  system  performance,  including  low- 
altitude  tracking  capability.  Netting  is  useful  in  several  respects  and 
at  several  levels  of  complexity.  First  of  all,  one  or  more  radar  sites 
t lacking  the  same  target  can  mix  (average)  their  vert ical-scan  error  sig- 
nals so  as  to  effectively  cancel  tracking-error  signals  due  to  reflections 
that  come  from  uncorrelated  patches  of  terrain,  while  the  direct  signals 
from  the  target  can  be  added  coherently.  This  is  a relatively  simple 
procedure  that  can  improve  low-angle  tracking,  especially  if  more  than 
two  sites  are  tracking  the  same  target,  and  is  useful  even  if  only  one  of 
the  sites  is  transmitting.  Multiple  transmissions  greatly  improve  the 
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tracking  capability,  however.  For  example,  consider  three 
sites  (on  separate  frequencies)  that  are  tracking  the  same 
exchange  information  about  return-pulse  signal  strength  in 


monopulse  radar 
target,  but 
the  elevation 


sum  and  difference  signals  derived  from  each  receiving  antenna 
each  frequency.  Let  the  parameters  be  defined  as  follows: 


a nd  for 


* Tmn  ~ TarK°t  reflection  coefficient  for  the  radar  beam  that 

is  reflected  from  transmitter  site  m to  receiver  site  n 

^mn  - Taiget  reflection  coefficient  for  the  radar  beam  from 

transmitter  m reflected  in  the  direction  of  the  reflect- 
ing surface  for  receiver  site  n 

pn  = Effective  surface  reflection  coefficient  between  the 
target  and  receiver  site  n 

n = Range  difference  between  the  direct  and  reflected 
signals  from  the  target  to  receiver 

* ^m  ~ sif°  m transmitter  frequency  (rad  s S 


- Range  from  the  target  for  the  signal  from  transmitter 
m,  as  received  by  site  n 


Assume  that  sufficient  knowledge  and  coordination  is  interchanged  between 
radar  sites.  Then  the  target  altituce  can  be  estimated  by  forming  a com- 
posite vertical-tracking  error  signal  to  drive  all  the  antenna  mounts 
together— i.e. , form  an  equivalent  error  signal  by  simple  averaging  of 
the  various  combinations  of  returns. 


To  consider  the  enhancement  in  low-angle  tracking  capability 

due  to  this  very  simple  (and  no  doubt  non-optimal)  averaging  of  the  various 

radar  signals,  assume  that  the  antennas  are  all  tracking  the  target  in 

azimuth  and  range  but  are  centered  on  the  horizon  (this  often  happens  in 

practice).  It  is  desired  that  an  estimate  of  the  target  elevation  angle 

be  made  so  as  to  generate  elevation  tracking  signals.  A single  monopulse 

radar  (see  Figure  15)  will  have  two  signals,  U and  L , derived  from  .he 

mn  mn 
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FIGURE  15  BLOCK  DIAGRAM  OF  A CONVENTIONAL  MONOPULSE  TRACKING  RADAR  Source:  Skoln.k  (1970.  p.  21-131 


upper  and  lower  horns,  respectively.  Let  A0  be  the  true  target  elevation 
angle  above  the  horizon  (and  lienee  above  the  radar  antenna  boresite). 

Now  if  we  assume  that  the  antenna  patterns  are  linear  functions  of  ele- 
vation angle  over  A0 , then  the  direct  return  received  by  the  upper  horn 
i s 
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where  (1  + A0)  is  the  antenna  pattern  magnitude 
phase  relative  to  boresight.  Thus  the  response 
to  both  the  direct  and  reflected  signals  can  be 
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and  e is  the  pattern 
of  the  upper  and  lower  horns 
written  as 
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Thus,  if  0 = keA0/2  and  1 - k , and  the  subscripts  m,  n are  understood, 
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If  ke/2  » 1 and  0 » A0  (and  since  sin  0 » A0  cos  0), 
we  have  a phase-comparison  monopulse  system.  The  estimate  of  A®  made 
by  this  system  is  ^6  given  by 
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If  ke/2  « 1,  then  an  amplitude-ccmparison  menupulse  system  results, 
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In  either  case,  if  a strong  reflection  occurs,  then  t|i  p 0 and  thus  0 st  1 1 
As  i varies  over  0 to  2tt,  A0  changes  from  0 to  infinity.  This  is  the 
essence  of  the  low-angle  tracking  problem. 


The  netting  process,  a simple  average  of  the  sum  and  dif- 
ference signals,  results  in  the  estimate  of  the  elevation  angle  A0  as 
given  by 
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In  the  case  of  a single  radar  system  (M  = N = 1),  it  is 
common  to  find  *P  « ^ and  hence  0 « 1.  Therefore,  as  71  changes,  A0  varies 
irom  near  zero  for  1 = 0 to  a very  large  value  when  7]  = jr  and  the  denomi- 
nator in  the  equation  above  approaches  zero.  So  for  a single  radar,  low- 
angle  tracking  is  often  impossible. 

Now  if  three  radar  sites  are  netted,  the  0 terms  add 

Eill  mn 

nill  Pmne  n approaches  zero.  The  result  of  summing 

the  0 terms  is  to  significantly  reduce  the  variance  of  A0  as  more  0 terms 

are  summed.  The  result  for  three  sites  and  hence  nine  paths  is  that  the 

variance  of  A0  is  reduced  from  something  like  A0  to  a few  percent  of  A0. 

This  will  produce  an  impressive  enhancement  of  low-angle  tracking  perfor- 
mance . 


Accurate  estimates  of  the  enhancement  could  be  made  if 

probability  density  functions  could  be  determined  for  p the  terrain 

n 

small-angle  reflection  coefficient.  As  mentioned  in  Section  II  above, 
sc.  ttering  theory  and  experimental  measurements  are  generally  lacking 
at  the  low  scattering  angles  of  interest.  However,  Rice  (1951)  has 
pointed  out  that  the  reflection  coefficients  (p)  measured  during  the 
survey  for  the  transcontinental  microwave  relay  system  could  be  repre- 
sented by  a Rayleigh  distribution  having  a median  value  of  about  0.28. 

7 ho  grazing  angles  at  the  specular  reflection  point  were  generally  less 
than  one-half  degree.  An  optimal  scheme  for  netting  radar  sites  should 
produce  some  improvement  over  the  simple  scheme  presented  here. 

If  each  radar  site  employs  one  or  more  of  the  other  low- 
altitude  improvement  techniques  discussed  in  this  report,  even  further 
improvements  can  result  from  netting  radar  sites  together.  Moving-target- 
indicator  (MTI)  clutter  rejection  is  not  deg- aded  by  netting  and  can  in 
fact  be  similarly  enhanced  if  desired. 
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, Passive  Radar  Receiver  Sites 


Some  of  'ho  advantages  of  nultistatie  operation  have  been 
covered  above.  The  remaining  consideration  of  multistatie  operation  is 
best  explained  by  examining  the  placement  of  a very  simple,  passive, 
nontracking  site  well  ahead  of  the  active  tracking  radar.  Such  a scheme 
is  illustrated  in  Figure  16,  Simple  time  of  arrival  (TOA)  of  the  direct 
return  from  the  target  and  the  target  return  as  received  by  an  omnidirec- 
tional and  pas  ' -e  receiver  can  then  be  used  to  determine  altitude  with- 
out the  need  to  derive  vertical  tracking-error  signals  from  the  active 
radar-site  signals.  Doppler  filtering  can  be  used  to  remove  clutter  in 
the  usual  manner.  Several  such  passive  receivers  might  be  used  around  a 
single  active  site  in  order  to  obtain  extensive  area  coverage  of  low- 
altitude  targets.  The  passive  sites  could  be  very  inexpensive  and  com- 
pletely unobservable  to  incoming  aircraft. 


TARGET 


FIGURE  16  RADAR  SYSTEM  USING  AN  ADDITIONAL  ANTENNA  (BISTATIC).  The  range 
and  azimuth  of  the  target  are  known  sufficiently  well.  The  elevation  is  determined 
by  the  time  of  arrival  of  the  radar  pulse  at  the  passive  site.  Several  such  sites 
could  be  used  to  obtain  a system  that  is  simple,  passive,  and  inexpensive. 


73 


The  effectiveness  of  such  a system  can  be  derived  by  con- 
sidering the  complexity  of  the  calculations  needed  to  determine  low-level 
target  altitude  and  by  dctci’mining  * he  area  coverage  of  a single  passive 
site. 


The  target  altitude  h^  is  easily  derived  from  the  following 

known  quantities:  the  target  range  (R  ),  the  passive  site  range  (R  ), 

T s 

the  difference  ir  azimuth  angles  to  the  target  and  to  the  passive  site 
(<3)  , and  the  range  between  the  target  and  the  passive  receiver  site  (r). 
Figure  17  illustrates  the  geometry  of  the  situation  for  which  the  fol- 
lowing equations  apply: 
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FIGURE  17  GEOMETRY  OF  THE  PASSIVE-RECEIVER  SYSTEM 
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(!. 


Frequency  Agility 


The  ability  to  dynamically  change  tie  radar  frequency  of 
operation  can  be  a great  aid  in  low-angle  tracking.  Barton  (1971)  gives 
an  excellent  detailed  discussion  of  this  technique  both  for  frequency 
hopping  and  spread  spectrum  when  signal  averaging  is  used,  and  only  a 
brief  additional  discussion  of  another  method  will  be  presented  here. 

The  usual  low-angle  tracking  problem  is  encountered  when  the  difference 
in  range  between  the  target  and  its  image  is  so  small  that  Doppler  filter- 
ing and  range  gating  cannot  affect  target  separation.  If  a frequency 
change  Af  in  the  transmitter  frequency  is  made,  then  the  relative  phase 
angle  between  the  target  and  its  image  will  change  by  2irAf6/c  radians 
where  6/c  is  the  range  difference  in  time  (on  the  order  of  0.1  ns).  Thus, 
an  extremely  wideband  radar  (/Af  > 1 GHz)  could  separate  the  target  and 
image  signals  by  appropriate  frequency  changes,  depending  on  the  received- 
si  gna  1 -to-noi  so  ratio.  These  dynamic  frequency-shifting  techniques  are 
closely  related  to  the  spread- spectrum  techniques,  but  may  be  considerably 
easier  to  implement  for  equivalent  results,  especially  if  only  accurate 
(racking,  as  opposed  to  area  scanning,  is  desired. 

One  method  to  accomplish  this  would  be  to  recursively 
estimate  the  target  range,  altitude,  and  ref lected-signa 1 phase,  and  from 
this  estimate,  the  altitude  difference  between  the  target  and  image.  Then 
the  next  probing  frequency  is  selected  to  change  the  relative  phase  so  as 
to  alternately  maximize  and  minimize  the  antenna  signal.  Altitude  infor- 
mation is  thus  not  derived  from  the  antenna  position  (which  centers  on 
the  horizon)  but  from  the  difference  frequency,  the  range  measurement, 
and  the  known  geometry  of  the  site. 
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6 ■ Multiple -Antenna  Arrays 

Multiple  antennas  provide  multiple  sources  of  information 
about  the  incoming  radar  echo  and  hence  may  be  used  to  implement  the 
netting  and  generalized  mot opulse  schemes  discussed  in  Sections  IJI-B-2-a 
and  III-B-2-b  above.  Multiple  antennas  may  be  used  in  a way  very  similar 
to  the  use  of  frequency  agility,  which  is  examined  in  Section  III-B-2-d 
above.  They  may  also  be  used  to  increase  the  effective  antenna  aperture 
as  shown  in  Section  IH-B-l-a  above.  Whether  or  not  additional  antennas 
can  provide  the  most  efficient  means  of  implementing  these  schemes  is  a 
question  of  engineering  tradeoffs  and  will  not  be  considered  here.  Our 
purpose  is  simply  to  show  possible  op  cions. 


In  the  case  of  netting  by  using  multiple  antennas  at  the 
same  site  it  is  necessary  to  show  that  the  antennas  used  sample  suffi- 
ciently different  reflected  echo  signals,  As  an  example,  consider 
antennas  displaced  vertically  at  the  same  site— i.e.,  having  different 
values  of  h^  m\ Figure  2.  For  the  netting  scheme  (Section  III-B-2-b) 
the  estimate  (Am)  of  the  true  elevation  angle  A0  is  given  by 
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For  the  scheme  to 
the  0mne  terms 


^n  ~ k^n»  where  k 


work,  the  values  of  T]n  must  be  different  enough  that 
add  incoherently  and  thus  tend  to  cancel  out.  Now 
= 2jtA  , and  6n  is  the  path  difference  between  the  direct 


t h 

signal  and  the  reflected  signal  to  the  n receiving  antenna.  The  question 

then  is,  does  «n  vary  sufficiently  rapidly  with  changes  in  receiving- 

antenna  height  (h^)?  The  approximate  expression  for  6 is  given  by  Kerr 

et  al.  (1951,  p,  415)  as  6 using  the  geometry  of  Figure  6.  The 

phase  difference  between  direct  and  reflected  echoes  is  T|  = 2^6 A,  so  the 

rate  of  change  of  T]  with  respect  to  receiver  antenna  height  (h  ) is 

r 

d~l/dh  = 4-th  /GX  . Considering  the  example  case  h = 100  m,  G = 10  km 
r t t ’ 

and  X = 3 cm,  we  find  dTj/dh^  = 4.2  rad  phase  change  per  meter  change  iu 
h^.  For  this  typical  case  we  see  that  indeed  T)  does  change  sufficiently 
rapidly  in  that  a pair  of  2-m-diameter  dishes  placed  one  right  above  the 
other  would  experience  more  than  an  8-rad  phase  change  between  them. 

The  frequency-agility  technique  (Section  III-B-2-a)  sepa- 
rates the  reflected  radar  echo  from  the  direct  radar  echo  by  noting  that 
the  phase  difference  between  the  direct  and  reflected  signals  changes  by 
an  amount  (2itAf6/c)  rad  when  the  frequency  is  changed  by  Af.  Holding  the 
frequency  constant  we  may  also  change  this  phase  difference  by  changing 
the  receiving-antenna  height  (h  ) . The  change  in  phase  difference  is 
given  by  (d^l/dlOAh^  = (4^h^/GX  )Ah^ . For  the  typical  case  mentioned 
above,  we  have  a change  in  phase  difference  of  4.2  rad  per  meter  change 
in  h^.  A disadvantage  of  changing  antenna  height  as  opposed  to  changing 
frequency  is  of  course  that  the  size  of  the  frequency  change  Af  can  be 
varied,  whereas  Ah^  would  have  only  a few  discrete  values  since  the  mul- 
tiple antennas  would  presumably  be  fixed. 

Multiple  antennas  could  very  well  be  used  in  the  formation 
of  multiple  antenna  beams,  as  is  done  with  array  antennas  (Skolnik,  1970, 
Chapter  11).  The  resui  irg  multiple  beams  would  then  be  used  with  a 
scattering  model  as  discussed  in  connection  with  the  generalized  signal- 
processing scheme  in  Section  III-B-2-a. 
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f 


Site  Location  and  Terrain  Modification 


f . 

The  elevation  error  induced  by  the  ground  reflection  of 
radar  echoes  generally  increases  as  the  reflection  coefficient,  p,  of  the 
surface  increases.  Hence,  one  method  of  improving  a radar  system's  low- 
angle  tracking  performance  is  to  reduce  p by  choosing  an  advantageous  site 
for  the  radar  or  physically  modifying  the  terrain  in  front  of  the  radar, 
or  both.  Many  radars  are  portable  and  can  be  moved  rather  easily  in  order 
to  improve  low-altitude  performance.  It  seems  likely  that  the  adverse 
effects  of  terrain  reflection  can  be  reduced  by  small  changes  in  radar 
site  location.  Once  a radar  site  is  chosen,  it  appears  likely  that  modi- 
fication of  a relatively  few  troublesome  spots  on  the  surface  can  be 
effective  in  further  reducing  p.  The  methods  that  can  be  employed  to 
reduce  P depend  strongly  on  the  scattering  mechanism  responsible  for  the 
reflection.  So  we  will  consider  specular,  quasispecular , and  diffuse 
scattering  mechanisms  separately.  Our  present  understanding  of  micro- 
wave  scattering  is  very  inadequate  in  many  situations  and  at  low  grazing 
angles  in  particular.  Therefore  a really  adequate  analysis  of  the  tech- 
niques suggested  here  awaits  new  advances  in  low-angle  scattering  theory 
and  experiments. 

The  specular  scattering  mechanism  is  responsible  for  the 
most  devastating  radar  errors  because  the  reflections  are  phase-coherent 

and  may  cancel  the  direct  echo  signal  completely.  The  specular  point  is 
generally  rather  close  to  a low-sited  radar,  making  modifications  of  the 
relevant  region  around  the  specular  point  fairly  convenient.  On  the  other 
hand,  for  low-angle  targets  the  first  Fresnel  zone,  from  which  most  of  the 
specular  reflection  comes,  is  rather  large.  For  example,  if  we  assume  a 
plane  geometry  as  in  Figure  6,  the  distance  tr  the  specular  point  X is 
only  about  0.5  km  from  the  radar  when  h^  = 5 m,  h^  = 100  m,  and  G = 10  km. 
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However  for  an  X-band  (\  = 3 cm)  radar  the  first  Fresnel  zone  is  an  ellipse 
extending  about  2.8  km  along  the  radar's  line  of  sight  and  about  30  m per- 
pendicular to  it.  As  noted  in  Section  II  the  specular  reflection  decreases 
markedly  as  the  surface  becomes  more  rough.  So  if  we  can  site  the  radar 
such  that  the  first  Fresnel-zone  region  about  the  specular  point  is  rela- 
tively rough,  specular  reflection  will  be  reduced  and  low-angle  performance 
aided.  If  the  site  is  fixed,  one  could  try  roughening  the  surface 

artificially  by  plowing  or  planting  appropriate  vegetation.  Evans  (1966) 

2 

reported  that  one-foot-high  grass  reduced  the  power  reflectivity  p of  a 

particular  site  by  a factor  of  16.  Another  idea  tried, out  by  Evans  (1966) 

was  to  erect  an  aluminum  screen  near  the  specular  point  to  shield  the 

specular-reflection  region  from  the  radar.  This  scheme  resulted  in  a 15-db 
2 

reduction  in  p for  the  most  favorable  elevation  angle  of  the  target.  As 
the  target  moves,  of  course,  the  specular  point  will  also  move.  Even  so, 
a single  screen  was  helpful  over  a range  of  elevation  angles.  Another 
option,  discussed  by  Smith  and  Melling  (1974),  is  to  have  the  ’’adar  sited 
high  enough  that  it  "looks  over"  the  specular- ref lection  region. 

The  quasispecular  scattering  mechanism  assumes  that  the 
reflecting  surface  is  broken  into  a number  of  small  reflecting  facets  or 
specular  points  (Kodis,  1966).  These  small  reflecting  regions  are  ran- 
domly tilted  with  respect  to  the  mean  vertical  direction  like  the  cracked 
surface  of  a boiled  egg.  Though  relatively  small,  the  facets  are  large 
compared  to  the  wavelength.  For  a given  geometry,  some  of  these  facets 
will  be  properly  oriented  for  specular  reflection,  and  it  is  these  facets 
that  provide  the  scattered  wave.  Usually  one  assumes  that  the  probability 
distribution  of  facet  tilts  is  Gaussian,  so  one  finds  fewer  properly 
aligned  facets  as  one  moves  away  from  the  specular  point.  The  properly 
aligned  facets  near  the  specular  point  provide  a coherent  surface  reflec- 
tion while  those  further  away  provide  an  incoherent  reflection  component. 
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As  the  average  tilt  of  the  facets  measured  by  the  unidirectional  rms 

slope  (tan  P ) increases,  the  incoherent  component  becomes  more  important, 
o 

As  with  specular  reflection  above,  siting  or  terrain  modification  to  make 
the  surface  near  the  specular  point  rough  will  reduce  the  coherent  compo- 
nent and  thus  aid  low-angle  radar  performance.  Once  the  region  around 
the  specular  point  is  sufficiently  rough,  one  expects  that  for  low-angle 
scattering  there  will  be  a relatively  small  number  of  properly  aligned 
facets  here  and  there.  If  these  troublesome  facets  could  be  located  and 
modified  either  by  earth  moving,  vegetation  cover,  or  possibly  conducting 
screens,  then  p would  be  reduced  still  further,  with  a resulting  improve- 
ment in  radar  performance.  Of  course  the  specular  point  and  troublesome 
facets  will  change  with  target  range  and  elevation;  but  as  with  the  shield- 
ing screens  mentioned  above,  it  is  expected  that  improvements  can  be 
effected  over  a range  of  target  parameters.  Virtually  no  closely  relevant 
experimental  data  have  yet  come  to  light,  and  this  is  just  what  is  needed 
to  determine  how  practical  these  schemes  are.  However,  Bullington  (1954) 
notes  that  the  horizontal  movement  of  an  antenna  by  only  100  ft  reduced 
the  low-angle  reflection  coefficient  p from  0,72  to  0.55.  Some  experi- 
mental suggestions  are  made  in  Section  IV. 

The  diffuse-scattering  mechanism  provides  an  incoherent 
reflection  from  virtually  the  entire  scattering  surface.  This  mechanism 
is  not  well  understood,  but  is  apparently  due  to  roughness  on  a scale 
small  compared  to  the  wavelength--e . g . , from  small  rocks  and  debris  on 
the  scattering  surface.  Diffuse  scattering  is  very  broad  and  not  strongly 
concentrated  in  the  specular  direction  as  are  the  specular  and  quasi  specu- 
lar scattering.  Hence,  the  (1/r^r^)  term  in  Eq.  (2)  is  important,  and 
diffuse  scattering  will  probably  arise  mainly  in  foreground  (near  the 
radar)  and  horizon  (near  the  target)  components.  So,  to  reduce  diffuse 
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scattering,  particular  attention  should  be  paid  to  the  foreground  com- 
ponent. Siting  the  radar  so  as  to  look  over"  the  foreground  component 
will  be  helpful.  Possibly  a layer  of  vegetation  over  the  surface  rocks 
and  debris  would  reduce  diffuse  scattering.  He  re  again,  relevant  theore- 
tical and  experimental  work  are  virtually  nonexistent  and  a critical 
evaluation  of  the  role  of  diffuse  scattering  in  low-argle  radar  perfor- 
mance awaits  new  developments. 

Once  one  finds  that  the  reduction  of  the  suiface  reflection 
coefficient  is  a practical  possibility  for  a given  situation,  then  one 
wants  to  know  just  how  much  low-angle  tracking  performance  is  improved 
for  a given  reduction  in  p.  This  calculation  is  rather  involved  and 
depends  mainly  on  the  following  parameters:  radar-target  geometry,  target 

elevation  relative  to  the  horizon  0^ , data  processing  used  at  the  radar 
(e.g.,  whether  or  not  data  smoothing  is  used),  radar  antenna  beam  patterns 
(sum  and  difference),  and  surface  scattering  model. 

Barton  (1974,  pp.  695-698)  has  looked  into  this  problem 
and  evaluated  the  expected  angular  error  fox'  various  values  of  the  afore- 
mentioned parameters.  For  a standard  monopulse  radar  he  finds  that  the 
transition  from  stable  to  unstable  tracking  occurs  when  the  target  eleva- 
tion angle  6^.  (see  Figure  11)  is  between  0.1  and  0.7  antenna  beamwidths 
(0e)’  dePendlnS  on  the  radar-tsrget  geometry,  the  surface  scattering  model 

and  the  radar  data-processing  £ lgorithm.  Whenever  p > 0.7  and  0 <070 

t # 0 * 

there  will  be  a strong  tendency  for  the  radar  to  track  a false  position 
near  the  horizon,  which  is  approximately  midway  between  the  real  target 
and  its  image  (see  Figure  1). 

It  is  difficult  to  generalize,  since  so  many  parameters 
are  involved.  However,  the  following  example  taken  from  Barton  (1974) 
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is  instructive  in  showing  the  sort  of  improvement  one  can  expect.  Con- 
sider a target  -t  10  km  range  being  tracked  by  a standard  monopulse  radar 
at  a height  of  5 m and  with  an  antenna  beamwidth  0 = 20  mrad.  The  rms 

elevation  error  (c^)  rises  as  the  target  position  approaches  the  horizon— 
i.e.  as  0t  approaches  zero.  Let  us  use  the  value  of  9 (let  us  call  it 
Q[)  at  which  the  error  ^ equals  0.2  0^  as  a figure  of  merit-the  smaller 
©t’  the  better  the  low-angle  performance.  Barton  analyzes  three  types  of 
scattering  surface:  a smooth  surface  with  rms  height  deviation  h =-  5 cm, 

a medium  surface  with  h = 25  cm,  and  a rough  surface  with  h = 1 m.  The 
values  of  9^  for  these  smooth,  medium,  and  rough  surfaces  are  15,  11,  and 
3 mrad,  respectively.  From  this  analysis  it  is  clear  that  substantial 
gams  in  low-angle  radar  performance  can  be  obtained  by  reducing  the  sur- 
face reflection  coefficient,  especially  near  the  specular  point. 

All  the  above  discussion  has  been  directed  toward  reducing 
P.  There  may  well  be  occasions  when  a particular  reflecting  facet  in  the 
horizon  region  could  be  used  along  the  lines  of  the  small  transponders 
discussed  in  Section  III-B-2-c  above. 


& • Multiple  Radar  Fences 

One  form  of  specialized  terrain  modification  is  the  implant- 
imr  of  radar  fences  to  block  the  return  signals  reflected  from  the  terrain. 
Some  benefits  could  result  from  the  proper  placement  and  adjustment  of 
multiple  fences  near  the  radar  site.  The  important  considerations  for  a 
single  fence  are  reviewed  by  Barton  (1974),  but  additional  benefits  from 
multiple  fences  were  not  considered.  In  a crude  way  of  thinking,  multiple 
fences  effectively  increase  the  antenna  aperture  so  that  the  main  antenna 
lobe  is  narrowed  near  the  horizon  enough  to  discriminate  between  the 
direct  target  and  the  reflected  returns. 
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Clutter  fences  are  commonly  used  to  reduce  backscatter 
from  the  nearby  ground.  In  certain  types  of  terrain  they  can  also  be 
used  to  block  be low-the-horizon  signals.  Such  a technique  has  been 
discussed  by  Hey  and  Parsons  (1955).  The  fences  can  be  placed  to  create 
an  artificial  horizon,  below  which  signals  will  be  strongly  attenuated 
as  shown  in  Figure  18,  The  major  disadvantage  of  the  multifence  horizon 
is  that  it  is  necessarily  flat,  and  in  uneven  terrain  it  may  be  possible 
for  a target  aircraft  to  fly  below  it  and  be  lost. 


ARTIFICIAL  HORIZON 
'VIAGE  OF  TARGET 

FIGURE  18  ARTIFICAL  HORIZON  CREATED  BY  MULTIPLE  RADAR  FENCES.  Thd  angle  0 
is  negative  below  the  artificial  horizon  and  positive  above.  The  image  of  the  target 
in  the  ground  is  suppressed  by  the  screens. 

In  order  to  get  a crude  estimate  of  the  effect  of  such  a 
row  of  fences,  we  shall  use  an  optical  approximation,  in  which  each 
fence  is  assumed  to  be  in  the  far  field  of  the  preceding  fence.  Let 
the  distance  from  the  target  to  the  first  fence  by  D,  and  the  distance 
between  the  first  fence  and  the  radar  dish  be  d.  Typically,  we  may  have 
D = 10  km  and  d = 100  m.  For  the  scattering  off  the  fence  nearest  the 
target,  we  can  assume  that  the  reflected  signal  from  the  target  is  a 
plane  wave  (since  D ->->  d).  The  scattered  intensity  is  given  by  the 
Fresnel  integrals  as  illustrated  in  Figure  19  (see  Born  and  Wolf,  1975, 
pp.  433-434). 
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FIGURE  19  SCATTERED  INTENSITY  AS  A FUNCTION  OF  ANGLE  6 ABOVE  (0  > 0)  AND 
BELOW  (0  < 0)  A SINGLE  KNIFE  EDGE  LOCATED  AT  6 = 0 

For  the  sake  of  illustration,  take  \ = 3 cm,  and  let  the 
diameter  of  the  dish  be  A = 2 m.  "hen  the  characteristic  fence  attenua- 
tion angle  9 zz/\/d  = 1°  is  comparable  to  the  main-lobe  width  of  the 
o 

dish,  9 = (1.22/)/A  = 1°.  The  main  lobe  will  be  reduced  slightly  on 

d 2 
the  ground  side,  and  the  lower  sidelobes  will  be  replaced  by  a 1/9 

monotonic  fall-off.  For  longei  wavelengths,  the  effect  of  the  fences 

on  the  main  lobe  becomes  more  significant. 

The  added  effect  of  the  additional  fences  is  difficult  to 

estimate,  but  a crude  guess  would  be  that  the  below-the-horizon  signal 
l N 

is  reduced  by  (4)  , where  N is  the  number  of  fences  in  addition  to  the 
outermost  fence. 

Although  the  below-the-horizon  signal  is  reduced  the 
fences  may  introduce  a new  error  for  above-the-horizon  signals  due  to 
the  oscillations  of  the  intensity  of  refracted  signal.  This  error  can 
be  reduced  in  several  ways: 

(1)  The  multiple  fences  can  be  spaced  in  such  a 
way  that  the  oscillations  from  the  separate 
fences  tend  to  cancel. 
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(2)  The  top  of  tho  fence  can  be  saw-toothed;  this 
should  have  the  effect  of  "apodizing"  the 
aperture  and  eliminating  the  high-spatial- 
frequency  oscillations. 

(3)  The  analysis  system  can  be  preprogrammed  to 
compensate  for  the  well  known  effect  that 
che  fences  will  have  on  any  particular 
radar  system. 

A more  precise  analysis  of  a multifence  system  will  require 
either  a more  detailed  mathematical  calculation  or  a computer  simulation. 

A final  evaluation  will  not  be  possible  until  such  work  is  completed. 

h.  Use  of  Azimuth  Information 

The  scattering  of  a radar  return  by  a rough  surface  occurs 
in  both  the  elevation  and  azimuth  directions.  Rough  terrain  also  has 
appreciable  correlation  lengths  in  both  directions,  so  that  the  scattering 
in  azimuth  should  exhibit  some  correlation  with  scattering  in  elevation 
(see  Figure  20).  This  phenomenon  cannot  be  used  alone  to  improve  low- 
altitude  tracking;  but  if  the  generalized  monopulse  technique  discussed 
in  Section  III-B-2-a  above  is  employed,  then  nonspecular  reflections  can 
become  the  major  difficulty  and  information  derived  from  the  azimuth 
tracking  circuits  can  be  used  to  estimate  and  hence  compensate  for  the 
diffuse  reflections.  A good  theory  of  radar  scattering  processes  would 
be  very  useful  in  quantifying  the  utility  of  these  signals. 

i . System  Optimization  Using  State-Space  Modeling 

All  the  techniques  discussed  in  this  section  have  uncer- 
tainties, and  difficulties  develop  for  weak  returns  and  extremely  low 
elevation  angles.  Thus  it  is  easy  to  argue  that  an  optimum  combination 
of  several  of  the  techniques  can  overcome  (or  at  least  improve)  some  of 
the  limitations  of  any  given  technique.  This  could  be  accomplished  in 
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practice  by  formulating  the  known  range  of  target  parameters,  the  radar 
characteristics,  etc.  into  a state-space  model  suitable  for  optimal  re- 
cursive estimation  of  the  target  altitude  as  a function  of  all  received 
signal  information.  The  block  diagram  of  such  a model  is  shown  in  Figure 
21.  This  would  be  a difficult  task,  both  in  the  initial  formulation  for 
a given  radar  configuration  and  in  the  real-time  calculations  required 
to  perform  the  Kalman  filtering,  but  it  could  have  a very  significant 
impact  on  radar  performance.  [For  reference,  Kalman  filter  is  discussed 
by  Schwartz  and  Shaw  (1975)].  If  an  existing  radar  network  were  to  be 
upgraded  in  capability,  this  additional  signal  processing  might  be  a 
very  attractive  alternative  to  other  major  changes  in  the  radar  system. 


FIGURE  20  CORRELATION  BETWEEN  SCATTERING  IN  AZIMUTH  (top  view)  AND  SCATTERING 
IN  ELEVATION  (side  view) 
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FIGURE  21 


GENERAL  MODEL  OF  RADAR  SYSTEM  SUITABLE  FOR  OPTIMAL  RECURSIVE 
ESTIMATION  OF  TARGET  PARAMETERS  SUCH  AS  HEIGHT,  RANGE,  AND  AZIMUTH 


IV  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Scattering  Theory 

1 • Present  Deficiencies  in  Scattering  Theory  Relevant  to 
Low-Angle  Radar  Tracking 

If  scattering  laws  for  a particular  range  of  parameters  are 
well  understood,  one  can  predict  with  confidence  how  the  direct  and  re- 
flected radar  echoes  will  differ,  and  algorithms  to  do  the  discrimination 
(see  Section  III)  can  be  evaluated.  Und  rstanding  of  the  scattering 
laws  may  also  suggest  discrimination  algorithms.  At  present  there  are 
many  deficiencies  in  sc.'tering  theory  that  manifest  themselves  as  un- 
certainties in  the  predicted  low-angle  tracking  errors  for  a given  radar 
system.  The  1960s  saw  considerable  advances  in  the  theory  of  electro- 
magnetic wave  scattering  from  statistically  rough  surfaces.  This 
relatively  rapid  advance  slowed  around  1970,  by  which  time  the  quasi- 
specular  theory  had  been  developed  and  had  satisfactorily  explained  a 
good  deal  of  experimental  daU  . However,  important  deficiencies  in  this 
and  other  theories  remain,  as  noted  in  the  following  paragraphs. 

Quasispecular  theory  breaks  down  at  low  grazing  angles  because 
shadowing  of  one  part  of  the  surface  by  another  and  multiple  scattering 
have  been  neglected.  Shadowing  theories  exist,  but  remain  largely  un- 
tested by  comparison  with  experiment.  A brief  comparison  of  a particular 
shadowing  theory  with  relevant  lunar  bistatie  radar  data  is  made  in 
Section  II  above.  Clearly  the  inclusion  of  shadowing  is  helpful,  but 
a comprehensive  comparison  remains  to  be  done.  No  treatment  of  multiple 
scattering,  which  is  particularly  relevant  at  low  grazing  angles,  has 
yet  come  to  the  authors'  attention, 
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Diffuse  scattering  is  thought  to  arise  mainly  from  surface 
roughness  on  scales  small  compared  to  the  wavelergth.  At  present  no 
really  good  statistical  or  geometric  theory  exists  for  the  diffuse 
component  for  any  grazing  angle.  The  semi-empirical  Lambertian  and 
Loirmel-Seeliger  scattering  laws  [Eqs.  (13)  and  (14)]  are  the  onlv  ones 
available.  To  date  no  comparison  of  these  laws  or  any  other  diffuse 
theory  appears  to  have  been  made  with  experimental  data  at  low  grazing 
angles. 

Scattering  from  vegetation  is  obviously  a difficult  pioblem, 
especially  at  low  grazing  angles.  According  to  Barrick  (1970),  the 
model  of  Peake  (1959a  and  1959b)  is  the  only  model  that  can  describe 
scattering  from  vegetation-covered  surfaces.  However,  Peake's  model 
does  not  agree  well  with  experimental  backscatter  measurements  at  low 
grazing  angles.  So,  for  the  case  of  special  interest  here,  we  must 
conclude  that  no  good  scattering  law  for  vegetation  exists. 

2.  Recommendations  for  Further  Work 

The  100-page  chapter  on  rough  surface  scattering  by  Barrick 
(1970)  in  the  ARPA-sponsored  Radar  Cross  Section  Handbook  was  very  use- 
ful in  the  present  study.  It  discussed  and  evaluated  existing  rough- 
surface  scattering  theory  up  to  1970  and  compared  theory  with  experimental 
results  for  backscatter.  It  is  recommended  that  a similar  effort  again 
be  sponsored  by  ARPA  to  include  work  since  1970,  with  an  emphasis  on  for- 
ward scatter  at  low  grazing  angles.  Such  an  effort  would  provide  an 
evaluated  collection  of  existing  theories  that  could  be  used  to  make 
interim  evaluations  of  radar  systems,  pending  advances  in  scattering 
theory.  It  would  also  provide  a good  jumping-off  point  for  new  theore- 
tical work.  The  substantial  research  results  in  the  Soviet  literature 
should  not  be  neglected.  A Russian  book  on  wave  scattering  from  a 
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statistically  rough  surface  by  Bass  and  Fuchs  (1972)  is  now  being  trans- 
lated and  edited  by  C.  B.  Vesecky  and  J.  F.  Vesecky  and  will  be  published 
by  Pergamon  Press  in  early  1976. 

The  deficiencies  in  scattering  theory  mentioned  above  still 
remain,  in  part  because  of  the  inherent  difficulty  of  the  low-angle 
scattering  problem  and  in  part  because  a great  deal  of  data — for  example, 
in  radar  astronomy — can  be  interpreted  successfully  without  worrying 
about  the  angles  near  grazing  incidence.  Therefore,  it  is  recommended 
that  ARPA  sponsor  theoretical  research  work  to  remedy  these  deficiencies 
at  low  grazing  angles.  There  are  probably  many  approximations  that  one 
can  exploit  at  near  grazing  incidence,  but  that  are  invalid  elsewhere. 
Therefore,  an  ARPA  program  could  well  benefit  from  emphasizing  the  de- 
velopment of  scattering  theory  that  may  be  valid  only  at  near  grazing 
incidence,  in  contrast  to  work  that  attempts  to  solve  scattering  prob- 
lems for  all  angles  of  incidence.  Quasispecular  scattering,  diffuse 
scattering,  and  scattering  from  vegetation  all  demand  attention  equally, 
though  it  might  be  possible  to  order  priorities  if  a specific  scattering 
terrain  were  known  to  be  of  special  interest.  Close  contact  should  be 
maintained  between  those  doing  theoretical  work  and  those  doing  the 
scattering  experiments  recommended  below. 

B.  Scattering  Experiments 

1 . Microwave  Experiments 

We  can  classify  microwave  experiments  into  two  broad  categories. 
First,  there  are  experiments  that  are  designed  to  be  critical  tests  of 
scattering  theory,  and  it  is  these  that  are  of  most  interest  here. 

Second,  there  are  experiments  designed  to  establish  empirical  scattering 
laws  for  specific  conditions.  There  are  a number  of  inherent  problems 
in  comparing  theory  with  available  experimental  data,  and  it  is  best  to 
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mention  these  problems  in  the  hope  that  future  experiments  will  seek  to 
avoid  them.  First,  experimental  results  are  often  quoted  in  terms  of 
the  Power  reflection  coefficient  (p2)  for  a particular  set  of  experimental 
conditions.  While  p2  is  a very  useful  quantity  in  practice,  it  is  not  an 
ideal  point  of  comparison  between  theory  and  experiment.  Since  p2  is  an 
integral  quantity  often  summing  up  contributions  from  large  portions  of 
the  scattering  surface,  agreement  between  theory  and  experiment  on  one 
or  several  values  of  o could  happen  by  chance  and  cannot  be  taken  as 
confirmation  of  a theory.  Experiments  should  note  the  variation  of  p 2 
over  a wide  range  of  experimental  parameters,  such  as  grazing  angle  and 
rms  surface  slope. 


Another  difficulty  arises  in  the  measurement  of  the  scattered 
radiation.  A full  specification  of  the  scattered  wave  requires  the 
measurement  of  the  general  Stokes  parameters  as  a function  of  direction 
and  time  of  arrival  as  well  as  Doppler  shift.  While  experimental  mea- 
surement of  a full  specification  is  a substantial  task,  it  may  be  neces- 
sary in  order  to  make  crucial  comparisons  of  theory  with  experiment. 

A genuine  test  of  theory  by  experiment  must  also  include  a 
sufficiently  accurate  characterization  of  the  scattering  surface.  For 
example,  quasispecular  theory  requires  measurement  of  the  unidirectional 
rms  surface  slope  on  a scale  larger  than  the  wavelength  k . Diffuse 
scattering  theory  would  probably  be  related  to  surface  variations  on 
scales  less  than  \.  So  surface  characteristics  need  to  be  measured  on 
scales  varying  from  much  smaller  to  much  larger  than  k . Often  the 
scattering  surface  during  an  experimental  run  is  so  poorly  known  in 
some  respects  that  the  resulting  data  cannot  be  compared  with  theory. 

It  is  best  to  choose  a scattering  surface  that  is  not  overly  complex, 

since  even  relatively  simple  surfaces  are  not  well  understood  at  low 
grazing  angles. 
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Designing  a good  experiment  that  allows  a crucial  comparison 
with  theory  will  not  be  an  easy  task.  The  following  is  a summary  of 
some  of  the  tasks  that  need  to  be  accomplished: 

(1)  Identify  key  assumptions  of  a theory  and  design 
experiments  to  find  out  when  these  assumptions 
break  down. 

(2)  Design  experimental  conditions  so  that  the  radiation 
scattered  by  different  scattering  mechanisms  may  be 
separated . 

(3)  Measure  a sufficiently  full  specification  of  the 
scattered  wave,  including  the  full  sot  of  general 
Stokes  parameters  if  necessary. 

(4)  Me  'sure  a sufficiently  full  set  of  scattering- 
surface  parameters  over  the  necessary  scale 
lengths,  which  may  be  both  much  longer  and 
much  shorter  than  the  wavelength. 

If  we  are  to  understand  rough  surface  scattering  at  low  graz- 
ng  angles,  it  will  be  necessary  to  conduct  experimental  work  to  help 
remedy  the  deficiencies  in  scattering  theory  noted  above.  It  is  there- 
foie  recommended  that  three  types  of  experiments  be  carried  out,  as 
described  below. 


a • Vegetation  Scattering  Experiments 

This  type  of  experiment  would  initially  be  directed 
toward  testing  the  Peake  (1959a  and  1959b)  geometric  scattering  theory. 
At  first,  vegetation  duplicating  Peake's  model  as  closely  as  possible 
should  be  used  as  a scattering  surface— i.e.  , vegetation  resembling  the 
dielectric  rods  randomly  distributed,  but  preferring  the  vertical,  and 
having  their  upper  ends  terminated  on  a horizontal  plane.  Microwave 
scattering  at  all  angles  of  incidence  would  be  explored,  but  with  a 
heavy  emphasis  on  grazing  incidence.  The  full  set  of  general  Stokes 
parameters  of  the  scattered  wave  would  be  measured  for  a variety  of 
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incident  polarizations.  The  experimental  objective  would  be  to  find  when 
and  why  the  model  breaks  down  by  varying  the  experimental  conditions,  one 
parameter  at  a time,  until  theory  and  experimental  measurements  diverge. 

b.  Quasispecular  Scattering  Experiments 

Quasispecular  theory  tends  to  break  down  lor  directions 
away  from  the  specular  direction  and  for  grazing  incidence.  The  experi- 
mental objective  would  again  be  to  find  just  when  and,  if  possible,  why 
the  model  breaks  down,  bv  varying  experimental  conditions.  The  scatter- 
ing surface  should  conform  as  closely  as  possible  to  the  gently  und'  fating 
surface  assumed  in  the  theory.  The  full  Stokes  parameters  of  the  scattered 
wave  should  be  measured  in  order  to  separate  out  the  quasispecularly 
scattered  radiation,  as  done  by  Tyler  and  Howard  (1973).  For  angles  near 
grazing  incidence,  the  shadowing  theory  of  Bass  and  Fuchs  (1972)  (see 
Section  II)  should  be  used  to  modify  the  quasispecular  cross  section  of 
Eq.  ( 1 1-4) . 

Another  interesting  quasispecular  scattering  experiment 
involves  a CW  transmitter  flown  at  low  altitude  over  the  scattering  sur- 
face. Such  an  experiment  using  the  Apollo  command  module  as  a transmitter 
and  an  Earth-based  receiver  observing  waves  scattered  off  the  lunar  sur- 
face is  described  in  Section  II.  A more  detailed  comparison  of  quasi- 
specular theory  including  shadowing  with  this  lunar  bistatic  data  (along 
the  lines  of  the  work  in  Section  II)  would  be  a useful  test  of  quasi- 
specular theory,  including  shadowing  near  grazing  incidence.  Another 
option  would  be  to  perform  a similar  experiment  with  a CW  transmitter 
flown  at  low  altitude  over  the  Earth's  surface.  Such  an  experiment  would 
in  principle  allow  one  to  find  the  number  and  location  of  the  "specular 
points"  responsible  for  the  quasispecularly  scattered  signal,  and  hence 
an  opportunity  to  try  out  the  terrain  modification  scheme  suggested  in 
Section  III. 
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Diffuse  Scattering  Experime n t s 


c. 

At  present  there  appears  to  be  no  really  useful  diffuse 
scattering  theory.  However,  the  semi-empirical  Lambertian  and  Lommel- 
Seeliger  scattering  laws  [Eqs.  (13)  and  (14)]  are  often  used  to  model 
diffuse  scattering  in  the  absence  of  better  alternatives.  It  would  be 
useful  to  know  just  how  well  these  laws  agree  with  experimental  measure- 
ments near  grazing  incidence  and  over  scattering  surfaces  of  varying 
roughness.  A well-conducted  experiment  paying  particular  attention  to 
surface  roughness  on  scales  small  compared  to  the  wavelength  could  well 
be  helpful  in  developing  a viable  diffuse-scattering  theory. 

2 . Laboratory  Experiments 

Although  one  naturally  thinks  of  microwave  scattering  experiments 
being  done  with  natural  terrain  in  the  field,  it  may  also  be  helpful  to 
conduct  experiments  with  carefully  prepared  model  surfaces  in  the  labora- 
tory. Near-field  and  edge  effects  would  complicate  laboratory  experiments, 
but  since  the  wavelengths  of  direct  concern  here  are  only  a few  centimeters 
at  most,  it  may  well  be  possible  to  overcome  the  difficulties.  In  addition, 
scaled  experiments  using  shorter  wavelengths  could  help  solve  these  prob- 
lems for  example,  by  using  millimeter-wave  sources  and  possibly  lasers. 
Laboratory  experiments  usually  provide  the  advantages  of  closer  control 
over  experimental  conditions.  Theories  based  on  laboratory  experiments 
would  still  have  to  be  tested  in  the  natural  environment. 

C . Radar  Systems 

1.  Data  Processing  Algorithms 

The  arguments  of  Section  III-B-2-a  demonstrate  that  many  of  the 
difficulties  that  monopulse  radars  encounter  during  low-angle  tracking 
are  simply  artifacts  of  the  particular  data-processing  algorithm  used. 
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With  appropriate  but  mode.st  changes  in  monopulse  radar  configuration 
and  data-processing  algorithms,  a radar  can  successfully  track  targets 
at  low  elevation  angles.  Practical  implementation  of  a first-order 
scheme  along  the  same  lines  as  suggested  in  Section  III-B-2-a  has  been 
accomplished  by  White  (1974)  with  very  encouraging  results,  especially 
over  water.  It  is  therefore  recommended  that  a program  bo  initiated  to 
accomplish  the  following: 

(1)  Study  a variety  of  radar  modifications  and 
data-processing  algorithms. 

(2)  Determine  the  sensitivity  of  the  above 
processing  algorithms  to  variations  in  the 
surface  scattering  model  and  other  radar 
system  parameters  such  as  signal-to-noise  ratio. 

■3)  Field  < ..  •cperimonts  that  demojtfstrate  the 

advantages  that  occur  with  modified  processing. 

The  theoretical  study  program,  (1)  and  (2),  would  be  at  the 
2-  to  3-man-year  level  for  one  year.  A comprehensive  theory  of  radar 
signal  processing  relative  to  multiple-feed  antennas,  antenna  arrays, 
and  radar  fences  and  involving  electromagnetic-wave  scattering  theory 
is  needed.  The  general  approach  could  begin  with  the  development  of 
optimal  algorithms  for  the  estimation  of  the  parameters  of  the  general 
model  presented  in  Figure  21.  Such  optimal  mathematical  programming 
algorithms  have  been  extensively  applied  in  control  theory,  but  very 
little  work  along  these  lines  appears  to  have  been  directed  toward  radar 
systems  and  none  at  all  toward  the  low-angle  tracking  problem.  An  im- 
portant input  to  the  theoretical  study  of  processing  algorithms  would  be 
the  results  of  the  studies  of  scattering  from  a rough  surface  at  grazing 
incidence  recommended  above. 
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2. 


mm-Wave  Radars 


Low-angle  tracking  performance  can  be  greatly  improved  in  a 
straight-forward  way  by  simply  reducing  the  radar  antenna  beamwidth.  Thus, 
the  interfering  multipath  signal  is  rejected  by  the  high  angular  resolution 
of  the  radar  antenna.  In  Sections  III-B-l-a  and  III-B-l-b,  two  methods  of 
achieving  sufficient  antenna  resolution  are  discussed.  The  first  is 
simply  to  use  a larger  antenna  aperture.  The  second  is  to  use  a higher 
operating  frequency.  The  numerous  advantages  of  small  antennas,  e.g., 
dishes  of  one  meter  and  smaller  diameter,  tend  to  make  one  favor  the 
latter  option.  With  regard  to  this  option  it  is  concluded  that  operating 
radars  at  higher  frequencies  (K  band  and  mm  wavelengths)  can  greatly  im- 
prove their  low-angle  tracking  performance  by  allowing  one  to  obtain  very 
narrow  antenna  beamwidths  with  relatively  small  antennas— the  higher  the 
frequency,  the  greater  the  improvement.  However,  as  one  operates  at  pro- 
gressively higher  frequencies,  progressively  higher  atmospheric  attenua- 
tion and  receiver  noise  levels  require  that  the  peak  transmitter  power 
(pt>  be  correspondingly  increased.  While  the  required  peak  transmitter 
power  at  10  mm  (K  band)  is  within  current  technological  capabilities, 
pt  rises  well  above  a megawatt  at  1,35  mm  for  a target  at  20  km  range. 

Since  current  transmitter  peak  powers  at  around  1 mm  fall  in  the  1-10 
watt  range,  substantial  increases  in  transmitter  peak  power  capability 
at  mm  wavelengths  are  required  in  order  to  exploit  the  superior  low-angle 
tracking  performance  at  these  wavelengths.  Improvements  in  receiver 
system  noise  temperature  over  the  estimates  used  in  Section  III-B-l-b 
could  also  be  helpful,  though  such  efforts  would  only  reduce  the  re- 
quired value  of  by  about  a factor  of  4.  It  should  also  be  noted 
that  at  short  wavelengths  (like  1.35  min)  antenna  surface  tolerances 
become  very  tight  indeed  (~0.05  mm)  if  satisfactory  performance  is 
to  be  achieved. 
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The  current  ARPA  program  in  high  power  mm  wave  sources  may  pro 
vide  the  necessary  high  power  levels  at  mm  wavelength.  When  mm  wave 
sources  in  the  megawatt  range  become  feasible,  research  and  development 
of  a mm  wave  low-angle  tracking  radar  should  be  considered  so  that  the 
superior  performance  at  mm  wavelengths  can  be  speedily  exploited.  As 
the  characteristics  of  high  power  mm  wave  sources  become  better  known, 
improvements  in  mm  wave  receiver  technology  should  be  studied  to  see  if 
funding  in  this  area  would  be  cost  effective  in  light  of  the  possible 
benefits  to  mm  radar  systems. 
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Appendix 

COMPUTER  PROGRAM  FOR  CALCULATING  RADAR  PARAMETERS 


BEGIN  COMMENT : 7 EST  LOW  ALTITUDE  RADAR  DI FFTIME . BY  DESPAIN  DEC  74 
FILE  IN(KI ND=REMOTE) ,OUT(KI ND=REMOTE) ; LABEL  HELL; 

FORMAT  DOUBLELINE 


FORMAT  SINGLELINE 
(" 


/-") ; 

REAL  R, I , A , RS , RM , TS , TM , AS , AM , FREQ , RMS , GAMMA  THETA , APERTURE , C , PI ; 

REAL  PROCEDURE  DI FFTIME (RANGE , TARGETHEIGHT , ANTENNAHEIGHT) ; 

REAL  RANGE, TARGETHEIGHT .ANTENNAHEIGHT;  BEGIN 
COMMENT : CALCULATION  OF  THE  TIME  DELAY  BETWEEN  THE  DIRECT 
RETURN  SIGNAL  AND  THE  REFLECTED  SIGNAL  FOR  LOW 

ALTITUDE  RADAR  TRACKING.  BY  DESPAIN  DEC  74 ; 

REAL  R , T , A , R2 , C , DELTA ; C : =299800000 ; A : =ANTENNAHEIGHT ; 
r. =TARGETHEIGHr;R2:=RANGE*RANGE; IF  (T+A)=0  THEN  DELTA : =0 
ELSE  DELTA: =SORT(R2+(A+T)**2)-SORT(R2+(A-T)**2) ; 

DI FFTIME : =DELTA/C; END  DI FFTIME; 

REAL  PROCEDURE  RHO (GAMMA , BETA) ; REAL  GAMMA , BETA ; BEG  IN 

COMMENT: J.  VESECKEY’S  BISTATIC  SCATTERING  FORMULA. BY  DESPAIN  74.; 

REAL  Q , RHOPRL , RHONML , E , P I , A ; 

REAL  PROCEDURE  LAMBDA(A) ; REAL  A ; BEGIN  COMMENT:  LAMBDA  FUNCTION; 
LAMBDA : =( 1 2*A)  * ( SORT ( 2/PI ) *EXP ( -A**2/2) -A*ERFC ( A/SORT( 2) ) ) ; END ■ 
PI : =0. 14159265359; 

E : =3 . 0 ; COMMENT : E IS  THE  SURFACE  DIELECTRIC  CONSTANT; 

IF  BETA  NEQ  0 THEN  BEGIN 
A: =TAN (GAMMA) / ( SORT (2) * TAN (BETA) ) ; 

Q : =( 1 /( 1+2+LAMBDA (A) ) ) ; END  ELSE  Q:=l; 

RHOPRL: =( (E*SIN(GAMMA)-SORT(E-COS(GAMMA)**2) ) / 

(E*S I N (GAMMA ) +SORT(E -COS (GAMMA) * * 2) ) )**2; 

RHONML: =( (SI N(GAMMA) -S0RT(E-C0S(GAMMA)**2) ) / 

(SI N (GAMMA ) +S0RT ( E -COS (GAMMA ) '* * 2 ) ) )**2; 

RHO: =. 5* (RHOPRL+RHONML) *Q ; 

END  RHO; 
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PI : =3. 141 o9265359;C: =299800000. ; 

Will  TE  (OUT,  <"ENTER  RADAR  CENTER  FREQUENCY (MHZ)  ">)  ; 

write  (out,  Center  rms  surface  slope">); 

WRITE(OUT,<-"ENTER  TARGET  HEIGHT  STEPS  (METERS  )">)  ; 

W R I TE  (OUT  , •'•"ENTER  MAXIMUM  TARGET  HEIGHT  (METERS)  ” > ); 
WRITE  (OUT,  -'"ENTER  ANTENNA  HEIGHT  STEPS  (METERS)  "»  ; 
WRITE  (OUT  , -/'ENTER  MAXIMUM  ANTENNA  HEIGHT  (METERS)  ">)  ; 
WRITE  (OUT  ENTER  RANGE  STEPS  (KILOMETERS)  ">)  ; 

WRITE (OUT,  ^"ENTER  MAXIMUM  RANGE (KILOMETERS )">)  ; 

F REQ : =FREQ* 1 . OftjOG ; RS : =1 000*  RS ; RM : =1 000*RM ; 


READ( IN , / , FREQ) ; 
READ (IN,/, RMS ) ; 
READ( I N , / , TS ) ; 
READ( IN, /,TM) ; 
READ( IN , / , AS) ; 
REA1)(  I N , / , AM)  ; 
READ( I N , / , RS ) ; 
READ ( I N , / , RM ) ; 


WRITE  (OUT  [.SPACE  4]); 

WRITE (OUT, -/’TESTING  THE  DIFFTIME  AND  THE  SCATTERING  LAW">) ; 

FOR  T : =TS  STEP  TS  UNTIL  TM  DO  BEGIN 
FOR  A : -AS  STEP  AS  UNTIL  AM  DO  BEGIN 
WR 1 TE ( OUT , DOU  BLEL I NE ) ; WR I TE (OUT [ SPACE  4 ] ) ; 

WR I TE (OUT , DOUBLELI NE ) ; WR I TE (OUT , 

<X10,"  LOW  ALTITUDE  RADAR  SCATTERING  PARAMETERS  "»  ; 

WR ITE (OUT , S I NGLEL I NE ) : WR I TE (OUT , < X 2, 

"TARGET  HEIGHT  =”,R6.2,”  METERS  , ANTENNA  HEIGHT  =”R6.2,”  METERS " 
>i T , A)  ; 

WHITE  (OUT,  SINGLELINE)  ; WRITE  (OUT  , <• 

RANGE  TIME  DELAY  SCATTER  RANGE  ANGLE  AMPLITUDE  APERTURE/ 

/»»  x 

>) ; 

WRITE  (OUT , <r 

(METERS)  (NANOSECONDS)  (METERS)  (RADIANS)  (RATIO)  (METERS)/ 

/"» ; 

WRITE (OUT, SINGLELINE) ; 

FYOR  R : =R 5 STEP  RS  UNTIL  RM  DO  BEGIN 
IF  R=0  THEN  GAMMA :=THETA:=1  ELSE  BEGIN 

GAMMA: =ARCTAN((A+T)/R) ; THETA : =ARCTAN( (T-A) / R) +GAMMA ; END : 

APERTURE: =C/( THETA* FREQ) ; 

COMMENT: STEGEN  FORMULA . REF: JASTK , "ANTENNA  HANDBOOK” , P2-24 ; 

WRITE  (OUTERS.  2 , R1 3 . 6 , R1  2 . 2 ,X2  , H9 . 6 , X2 , R9 . 6 , R1 1 .4  >, 

R , DI FFTIME(It  ,T , A)*0+9  R*A/(A+T) , THETA , RHO (GAMMA , RMS ) .APERTURE) ; 

END; END; END; 

HELL : END  TES TD 1 F FT I ME . 
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